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ABSTRACT 


This  report  is  concerned  with  the  development  and  presentation  of 
orthogonal  main-effect  plans.  These  plans  permit  uncorreluted 
estimates  of  all  main  effects  of  both  symmetrical  and  asymmetrical 
factorial  ejcperiments  with  a  minimum  number  of  trials. 

Chapters  II  and  III  outline  some  background  material  on  fitting  linear 
models  and  factorial  experiments  which  the  user  of  this  report  may 
find  informative.  These  tv  .  chapters  give  a  short  review  of  existing 
knowledge  of  factorial  experiments  and  methods  of  analysing  them. 

Chapter  IV  gives  an  account  of  the  development  of  orthogonal  main- 
effect  plans  for  symmetrical  and  asymmetrical  factorial  experiments. 

The  plans  for  asymmetrical  experiments  are  based  on  the  proposition 
that  if  the  levels  of  a  factor  occur  with  the  levels  of  another  factor  with 
proportional  frequencies  then  the  two  factors  are  orthogonal.  The 
possibilities  of  blocking  these  plans,  the  efficiencies  of  the  estimates, 
the  randomization  procedure  and  the  method  of  analysis  are  discussed. 

The  report  concludes  with  a  catalogue  of  orthogonal  main-effect  plans. 
This  catalogue  consists  of  the  treatment  combinations  of  twenty-six  basic 
plans,  involving  factors  with  up  to  nine  levels  and  with  up  to  eignty-one 
trials,  from  which  all  orthogonal  main-effect  plans  which  can  be 
constructed  with  eighty-one  or  fewer  trials  may  be  deduced. 
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I.  INTRODUCTION 
A.  Preliminary  Remarks 

The  purpose  of  this  report  is  to  present  in  as  simple  form  as  possible 
a  catalogue  of  plans  by  which  the  effects  of  ?  number  of  controllable 
variables  can  be  explored. 

The  general  situation  is  that  there  are  response  or  resultant  variables 
or  outputs  which  are  thought  to  depend  on  controllable  variables  or  inputs, 
as,  for  instance,  the  response  of  a  chemical  process,  which  is  a  resultant 
variable,  depends  on  temperature  of  reaction,  pressure,  type  of  catalyst 
present,  flow  rate  of  ingredients  and  so  on.  The  situation  is  one  of  very 
general  occurrence  as  may  be  seen  from  the  following  examples  from 
widely  different  areas  of  human  investigation. 


Situation  Response 

1.  Performance  of  college 
freshmen  students  in 
mathematics 

2.  Performance  of 
astronauts 

3.  Conversion  of  one 
chemical  to  another 


Inputs 

Number  of  lectures,  method 
of  presentation,  number  of 
assignments 

Types  and  qualili-’s  of  foods, 
other  possible  environment 
stimuli,  amount  of  training 

Temperature,  pressure, 
feed  rates,  catalysts, 
contact  time 
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Situation 

Response 

Inputs 

4. 

Quality  of  an  electrical 

device 

Variables  in  method  of 

production  of  device,  such 

as  nature  of  alloy,  of 

resistances,  rate  of  cooling 

in  production  of  parts 

5. 

Growth  of  a  t-‘ological 

organism 

Amounts  and  types  of 

various  nutrients 

6. 

Yield  of  an  agricultural 

crop 

Rate  of  seeding,  spacing  of 

plants,  amounts  of 

fertilizers 

7, 

Psychological  status  of 

sick  individuals 

Amounts  of  drugs,  amount 

and  nature  of  psychoanalysis 

8. 

Treat. nent  of  an  illness 

Diet  factors,  drug  factors, 

amount  of  rest 

V. 

Degree  of  delinquency 

of  humans 

Social  and  economic 

measures 

The  JL  eason  for  making  a  list  like  the  above  which  could  be  extended 
indefinitely  is  to  show  the  range  of  situations  which  have  the  same 
essential  st^uc^•^^e.  From  the  point  v«ew  of  designing  an  investigation 
in  any  of  these  situations  the  problems  are  as  follows: 

(i)  Defining  in  operational  terms  the  resultant  or  response  variables 
of  interest.  In  the  above  examples  only  one  response  variable  is 
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given,  but  one  can  easily  imagine  many  others. 

(ii)  Defining  the  control  variables,  or  inpuis  which  should  be 

considered.  For  various  reasons  statisticians  referred  to  these 
inputs  as  factors  and  this  led  to  “he  term  factorial  experiments, 
■which  are  nothin  •  but  experiments  de.signcd  to  investigate 
several  factors  or  inputs. 

(iii)  Defining  the  variants  of  factors  to  be  considered,  as  for  instance 
temperature  at  150*^C,  Z00*^C.  Z50°C,  or  catalyst  as 
manganese  or  platinum  oxide.  It  is  fairly  standard  to  use  the 
-erm  "levels'-  for  these  variants.  In  the  uase  of  a  factor  like 
temperature  which  can  be  envisaged  rs  taking  any  value  in  a 
partic*ular  range,  that  is,  a  continuous  factor,  the  term  "le-vel"  is 
clearly  appropriate.  For  the  case  of  discrete  factor,  or  one  Lr. 
which  the  variants  cannot  be  represented  as  points  on  a  line,  the 
term  "levels"  is  not  as  appiopriate  because  one  variant  cannot  be 
said  to  be  at  a  higher  level  than  another.  But  the  use  of  the  term 
"levels"'  for  both  cases  does  not  appear  to  be  confusing  to 
scientists  and  technologists,  and  is  so  well  entrenched  in 
statistical  t<  rmin'>logy,  that  it  will  be  followed  here. 

(iv)  Specification  of  the  class  of  situations  to  which  the  inputs  are  to 
be  applied  and  about  which  conclusions  are  desired. 

{v)  Choice  of  combinations  of  the  inputs  to  be  tested. 

(vi)  Assignment  of  the  individual  combinations  to  the  members  of  the 
class  of  situations. 
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(v'ii)  Specification  of  methods  of  interpreting  the  resultant  data. 


All  the  above  are  shorthand  statements  of  problems  each  of  which  could 
be  given  extensive  consideration.  Some  of  the  problems  are  strictly  in 
the  province  of  the  experimenter.  The  statistician  per  se  does  not  know 
what  yield  variables  are  of  interest,  what  possible  inputs  should  be 
considered,  how  the  yield  variables  are  to  be  measured,  how  the  inputs 


are  to  be  controlled,  what  levels  of  the  inputs  are  to  be  considered,  and 
what  class  of  situations  is  of  interest  to  the  experimenter.  The 
statistician  can  sometimes  give  advice  on  these  matters,  based  on 


[■atiuiis  uf  repeatabllliy  of  ability  to  coiiti'ol  vaiTabies,  ot  pre¬ 


cision  and  sensitivity  of  possible  choices,  variability  c£  members  of  the 
class  of  situations  to  be  investigated,  and  precision  of  measurement  of 
yield  variables.  To  evaluate  whether  a  statistician  can  give  help  on  these 
matters,  he  will  ask  questions  of  the  experimenter  pointing  out  the 
consequences  of  various  conditions  and  choices  and  may  on  the  basis  of 
the  answers  make  suggestions  or  merely  content  himself  with  the  opinion 
that  the  experimenter  has  already  considered  these  aspects  adequately. 

It  is  when  we  turn  to  the  latter  problems  of  the  above  list  that 
particular  knowledge  of  statistical  design  of  experiments  comes  into  play. 
One  may  for  instance  envisage  a  situation  in  which  there  are  say  10  inputs 
or  factors,  each  of  which  could  be  examined  at  several  levels.  The 
naive  reaction  is  to  say:  "Try  all  possible  combinations",  but  when  one 
realizes  that  even  if  the  number  of  levels  chosen  for  each  factor  is  3 
the  total  number  of  possible  combinations  is  3^^  or  59,049,  one  sees 
that  thi*"  is  completely  impossible.  Also  it  may  not  be  necessary.  To 


see  this  we  need  to  consider  further  the  type  of  problem  being  attacked. 
That  problems  can  be  classified  is,  we  imagine,  self-evident,  but 
possibly  one  of  the  real  gains  from  statistical  thinking  is  the  existence 
of  a  classification,  which  can  of  course  be  only  rough. 

It  will  be  convenient  in  what  follows  to  use  sometimes  the  phrase 
■'factor  space".  This  is  a  short  term  for  the  totality  of  possible 
combinations  of  the  factors  or  inputs  which  is  considered  relevant.  If 
for  instance  one  wishes  to  investigate  temperature  between  150°  C  and 
300°C,  pressure  between  lOlbs/sq.in.  and  20  Ibs/sq.  in. ,  flowrate 
from  100  gallons  per  minute  to  500  gallons  per  minute  with  no 
restrictions  as  to  what  combinations  of  temperature,  pressure  and  flow 
rate  are  possible,  the  facto.r  space  can  be  represented  as  a  rectangular 
parallelepiped  in  3  dimensions  with  perpendicular  axes  representing  the 
three  factors.  The  interior  of  this  parallelepiped  then  contains  a 
representation  of  all  possible  combinations  of  factors. 

Even  when  the  ’xesponse  variables  or  output  variables,  and  the 
control  or  input  variables  have  been  defined,  problems  can  be  classified 
according  to  the  end  result  desired; 

{A)  The  aim  may  be  to  determine  the  combination  or  combinations  of 
input  variables  which  gives  maximum  response  or  minimum 
response.  This  is  entirely  obvious  in  the  case  of  a  chemical 
process  in  industry,  or  in  the  training  of  a  skilled  operator.  For 
brevity  this  problem  is  referred  to  as  the  problem  of  response 
optimization.  There  may  be  several  response  variables  which  are 
to  be  optimized  jointly  and  one  may  then  get  into  problems  of 
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linear  or  non-linear  programming  as  well. 


(B)  The  aim  may  be  to  determine  which  of  the  possible  factors,  which 
can  be  imagined  as  possibly  affecting  response,  do  indeed  have  a 
non-trivial  effect.  If  one  has  a  production  line  involving  many 
distinct  stages  which  is  producing  articles  which  are  not  acceptable, 
one  has  the  problem  of  determining  which  of  the  possible  variables, 
of  which  there  can  easily  be  20  or  more,  are  affecting  the  quality 

of  the  end  product.  T^is  problem  will  be  referred  to  as  the 
problem  of  screening  of  factors. 

(C)  The  aim  may  be  to  obtain  a  rough  idea  of  the  effects  of  factors 
applied  jointly  over  a  range  of  conditions.  This  problem  has  nut 
had  a  particular  name  associated  with  it,  and  for  lack  of  a  better 
term,  we  shall  call  it  factorial  evaluation,  that  is,  evaluation  of  the 
role  of  the  possible  factors. 

(D)  The  aim  may  be  to  obtain  a  continuous  functional  relationship  of  the 
response  (output)  to  the  factors  (inputs)  like,  for  instance, 


where  y  is  the  response,  and  Xj,  are  the  inputs.  This 
problem  will  be  referred  to  as  functional  evaluation. 

It  is  not  our  purpose  here  to  discuss  all  these  problems  in  detail,  but 
the  follo'.ving  remarks  indicate  to  some  extent  what  is  involved  and  what 
general  plan  of  investigation  should  be  followed. 

In  the  case  of  response  optimization,  it  i^  not  essential  to  obtain 


6 


much  of  an  idea  of  what  factors  are  relevant,  provided  one  has  included 
in  one’s  list  all  the  factors  which  are  controllable.  One  can  merely  do 
what  might  be  termed  local  experimentation  in  the  factor  space,  by 
experimenting  around  a  point  in  the  factor  space  which  is  thought  to  be 
the  best  guess  of  the  optimal  combination.  On  the  basis  of  this  experi¬ 
mentation  one  finds  the  direction  in  the  factor  space  along  which  it  seems 
best  to  p-  :  ceed  in  search  of  the  optimum.  The  formulation  of  a  strategy 
by  which  to  do  all  this  is  a  difficult  problem,  but  some  statisticians  have 
in  recent  years  put  forward  some  interesting  ideas.  It  is  of  course 
apparent  that  the  more  factors  considered  the  greater  the  complexity  of 
assessing  the  experimental  results.  If  there  should  be  two  or  more 
response  criteria,  the  problem  of  optimum  seeking  becomes  much  more 
complicated.  For  instance  one  might  wish  to  determine  the  combination 
of  inputs  which  gives  maximum  percentage  response  of  a  chemical 
process  subje-t  to  the  restriction  that  one  wishes  to  attain  a  purity 
greater  than  a  certaL  percentage,  with  as  low  a  use  of  catalyst  as 
possible.  Some  of  the  problems  which  occur  on  first  thought  may  not 
even  be  sufficiently  well-defined  to  enable  even  the  theoretical  search  for 
a  solution.  They  may  also  depend  on  obtaining  a  fairly  good  idea  of  the 
functional  structure  of  the  situation,  that  is,  the  mathematical  relation¬ 
ships  of  the  responses  or  outputs  to  the  inputs. 

In  the  case  of  the  experiment  for  screening  factors,  the  problem  is  more 
one  of  determining  factors  which  are  having  appreciable  effects  with  a  view 
to  more  precise  experimentation  directed  to  aims  j(A),  (C)  or  (D).  For 
instance,  suppose  that  a  certain  step  in  a  production  process  consists 
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of  heating  the  uncompleted  product  entering  that  stage  in  a  furnace  for 
two  hours  at  1S0°C.  One  might  ask  whether  it  matters  whether  the 
heating  is  don^  for  a  much  shorter  time  like  half  an  hour  or  a  longer 
time  likt  -  hours.  In  other  words  is  this  a  factor  which  merits  some 
detailed  e>.amination  or  can  one  assume  that  realistic  changes  in  the 
factor  level  are  going  to  have  an  inappreciable  effect  on  the  resultant 
product.  A  procedure  commonly  used  for  this  sort  of  investigation  is  to 
use  the  following  plan  in  wh'ch  there  are  6  factors  and  L  denotes  a  low 
level,  K  a  high  level  (or  in  the  case  of  a  discrete  factor,  two  interesting 
possibilities): 

Factor 


Trial 

1 

2 

3 

4 

5 

6 

1 

L 

L 

L 

L 

L 

L 

Z 

H 

L 

L 

L 

L 

L 

3 

L 

H 

L 

L 

L 

L 

4 

L, 

L 

H 

L 

L 

L 

5 

L 

L 

T 

H 

L 

L 

6 

L 

L 

L 

L 

H 

L 

7 

L 

L 

L 

L 

L 

H 

It  is  not  at  all  difficult  to  make  a  convincing  case  that  this  plan  is  very 
■lefficient.  The  problem  is  to  evaluate  6  factors  and  the  later  parts  of 
this  technical  report  exhibit  plans  which  can  be  shown  to  have  greatest 
efficiency  and  sensitivity  in  determining  whether  factors  merit  further 
study.  These  are  the  main-effect  plans  for  which  a  catalogue  is  given. 
The  development  and  cataloging  of  these  plans  were  the  main  object'ves 
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of  the  present  research. 

When  we  turn  to  what  we  have  termed  "factor  evaluation",  we  are 
interested  in  not  only  what  factors  have  effects  of  non-trivial  importance, 
but  whether  also  the  effect  of  one  factor  depends  on  the  status  of  other 
factors  at  which  this  effect  is  determined.  In  standaid  statistical  jargon 
the  question  is  "What  are  the  effects  and  interactions  of  the  factors?"  . 
For  this  sort  of  task,  the  gamut  of  factorial  experimentation  as  developed 
over  the  past  30  years  is  rel'vant. 

Finally  when  we  consider  the  evaluation  of  functional  relationships  we 

not  only  want  to  know  what  factors  and  interactions  are  present  but  we 

want  to  express  the  relationship  in  as  scientifically  meaningful  way  as 

possible  and  we  have  to  take  account  of  the  units  in  which  factor  levels 

are  measured,  and  have  to  search  for  underlying  variables  which  may  be 

composites  of  the  variables  on  which  we  choose  to  experiment.  For 

instance  we  may  experiment  on  a  variable  v  which  is,  say,  a  velocity, 

but  the  way  velocity  enters  into  the  determination  of  response  is  in 
1  ^2 

terms  of  (v  +  b)  ' 

Tliore  are  common  elements  to  all  these  aims  and  there  are  no  sharp 
divisions  among  them.  In  many  cases  finding  the  optimiun  is  the 
ultimate  aim,  but  screening  of  factors  and  looking  for  the  possible 
existence  of  interactions  is  undertaken  first.  Similarly  screening  of 
factors  and  evaluation  of  interactions  may  well  precede  the  search  for  a 
functional  relationship.  So  the  approach  to  a  problem  of  science  or 
technology  is  a  matter  of  judgment.  An  aim  of  the  theoretical  study  of 
design  of  experiments  is  to  construct  a  rationale  to  aid  the  reaching  of 
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such  a  judgment. 


B.  General  Background  ot  Material  Presented 

The  aim  of  the  research  underlying  t;‘;s  report  is  to  present  a 
catalogue  of  plans  which  v/ill  enable  the  experimenter  to  screen  factors. 
The  plans  enable  the  estimation  of  the  effect.s  of  all  the  factors  included.. 
Any  such  estimation  is  unbiased  if  there  are  no  Interactions.  If  there  are 
interactions  estimates  obtained  by  a  model  assuming  absence  of  inter¬ 
actions  will  be  deviate  from  their  true  values  by  other  than  ejcpeiimental 
error.  This  should  not  be  regarded  as  a  deiiciency  ot  the  rSa'i-,  because 
the  essence  of  research  is  the  obtaining  of  ideas  which  are  subjected  to 
confirmation.  To  demand  that  an  experiment  have  a  completely 
unambiguous  interpretation  is  realistic  only  if  the  experiment  will  not  be 
repeated,  that  is,  if  it  is  a  tenninal  one,  and  such  experiments  must  be 
rare.  No  decisions  in  research  are  irreversible,  and  knowledge 
possessed  at  a  particular  point  of  time  is  at  best  an  approximation  to  the 
truth  and  at  worst  completely  fallacious.  Questions  underlying  this  state¬ 
ment  can  easily  be  fotmulated,  and  one  may  question,  for  instance,  the 
risks  involved  in  any  plan  of  investigation. 

In  addition  to  the  non-terminal  nature  of  research  conclusions,  one 
niust  also  take  into  account  what  might  be  termed  the  economy  of  research. 
One  can  envisage  using,  at  a  particular  stage  of  an  experimental  investi¬ 
gation,  a  range  of  plan?  from  the  smallest  and  least-time-consuming  plan 
which  will  enable  one  to  get  some  ideas,  to  a  large  expensive  plan  which 
will  give  clear-cut  unambiguous  .answers.  With  the  former  there  is  the 
risk  of  reaching  erroneous  conclusions,  but  the  .advantage  of  getting  a 
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rough  picture  quickly.  With  the  latter  the  ris.k  of  reaching  erroneous 
conclusions  will  be  low,  but  the  chance  of  reaching  conclusions  which  are 
highly  uninteresting  may  be  quite  appreciable.  Also  if  the  experiment  is 
to  take,  say,  3  months  to  perform,  one  may  well  find  that  the  ideas  which 
led  to  its  being  planned  have  been  modified  by  experience  and  knowledge 
acquired  since  the  planning,  so  that  the  "big"  experiment  only  partially 
done  is  clearly  inappropriate  and  misdirected.  In  the  case  of  technologi¬ 
cal  experiments  in  industry  there  is  obviously  a  value  to  be  gained  from 
approximate  conclusions  obtained  quicldy.  Even  in  what  might  be  termed 
pure  research  of  no  conceivable  economic  or  social  value,  the  researcher 
will  be  concerned  about  the  utilization  of  his  own  time  and  energy.  It  is 
apparent  that  one  should  commit  oneself  to  a  large  experiment  which  is 
seeking  a  detailed  picture  only  after  one  has  identified  factors  or  inputs 
which  are  known  to  have  interesting  effects  and  interactions. 

The  catalogue  is  then  a  catalogue  of  experimental  plans  which  are 
likely  to  be  useful  in  ixploratory  research.  The  adieCtive  "exploratory" 
here  is  not  meant  to  imply  research  based  on  little  ki.owledge  but  research 
perhaos  in  an  area  which  is  highly  developed,  where  one  wishes  to  obtain 
a  quick  idea  of  which  factors  whould  be  investigated  mere  deeply  and 
wh  ich  factors  should  be  ignored.  There  are  of  course  risks  involved  in 
ignoring  a  factor  or  in  deciding  that  variation  in  levels  of  that  factor  is 
iiv/t  worth  including  in  the  investigation.  It  may  be  that  the  factor  has  an 
interesting  effect  in  only  a  small  range,  as,  for  example,  a  biological 
stimulus  such  as  an  estr-^gen.  For  example,  it  was  known  for  years 
that  stilbestrol  caused  some  species  of  animals  to  have  increased 
growth  rates,  but  it  was  found  that  with  doses  which  were  thought  to  have 
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any  possible  effect  the  side  effects  were  intolerable-  Later  it  was  found 
that  doses  which  were  small  relative  to  doses  previously  tried  had  the 
desired  effects  with  none  of  the  undesired  ones. 

There  is  some  further  insurance  of  uncertain  value  in  the  use  of  these 
plans,  which  arises  from  the  empirical  conclusion  that  there  are  not 
likely  to  be  sizeable  interactions  if  there  are  no  main  effects.  This  does 
emphasize  that  one  should,  by  one  way  or  another,  have  some  check  on 
the  magnitude  of  error  in  ne  situation  being  examined,  because  the 
determination  of  whether  there  are  effects  of  interesting  magnitude 
depends  on  two  things  (1)  whether  the  actual  numericed  magnitude  is 
interesting  and  (ii)  whether  the  actual  magnitude  is  sizeable,  say  of  the 
order  of  1-^  or  2  times  its  standard  error. 

The  catalogue  of  plans  enables  an  experimenter  to  discover  quickly 

what  pJans  are  avail  ible  for  his  particular  situation.  He  may  for  instance 

wish  to  look  at  two  factors  at  five  levels,  three  factors  at  four  levels, 

two  factors  at  three  levels  and  one  factor  at  two  levels.  In  the  technical 

language  common  to  the  area  of  the  design  of  experiments,  he  is  involved 
2  3  2 

in  a  5  x4  x3  x2  factorial  situation.  To  list  all  possible  plans  would 
be  an  impossible  task  and  we  have  confined  ourselves  to  plans  which 
require  no  more  than  81  observations.  The  plans  listed  are  orth-ogonal 
ones,  that  is,  they  enable  best  unbiased  estimates  of  effects  of  all 
factors  which  are  uncorrelated.  Even  to  set  out  all  the  possibilities  in 
this  case  would  be  tedious  but  sonie  condensation  of  the  listing  is 
accomplished  by  giving  an  index  with  instructions,  so  that  plans  can  be 
used  with  minor  modifications  for  other  situatior  s. 
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One  modification  of  standard  plans  which  is  always  possible  has  been 
little  used  in  the  past.  This  modification  was  used  in  the  construction  of 

the  plans  and  can  be  used  to  a  wider  extent.  If  we  have  a  situation  like  a 

.-2  3  2  8 

5  X  4  X  3  X  2  mentioned  above  we  can  use  a  plan  for  a  5  experiment 

and  replace  three  five-level  factors  by  four-level  factors,  two  five-level 

factors  by  three-level  factors,  and  one  five-level  factor  by  a  two-level 

factor.  In  the  last  case  one  would  set  up  the  following  correspondence: 

level  of  five-’ avel  factor  0  12  3  4 

level  of  two-level  factor  0  1110 

Thus  levels  0  and  4  of  the  five-level  factor  arp  replaced  by  the  0 
level  of  the  two-level  factor  and  levels  1,  2  and  3  of  the  five -level 
factor  by  the  1  level  of  the  two-level  factor.  If  one  really  wanted  to 
experiment  with  some  six-level  factors  one  could  collapse  a  seven-level 
factor  plan.  This  results  in  a  little  loss  of  statistical  efficiency,  bu*  not 
enough  to  worry  about.  At  least  it  seems  preferable,  to  the  present 
authors,  to  encounter  a  small  loss  in  efficiency  in  order  to  accommodate 
the  six-level  factor  rather  than  to  force  the  experimenter  to  delete  one  of 
the  levels  he  likes  or  otherwise  revamp  the  situation.  Of  course  there  is 
no  point  in  introducing  levels  merely  'cr  the  sake  of  doing  so,  and  the 
more  levels  that  are  included  for  a  particular  factor,  the  more  trials 
a  e  required. 


C.  Structure  of  the  Material  Presented 

The  analysis  of  the  orthogonal  main-effect  plans,  i.  e.  estimation  of 
parameters,  estimation  of  error,  tests  of  significance,  is  the  standard 
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one  based  on  the  method  of  least  squares  and  a  brief  account  of  the 
features  of  this  method  is  given  in  Chapter  II. 

The  basis  for  most  of  the  plans  is  the  concept  of  factorial  experi¬ 
mentation  and  the  elementary  ideas  of  this  topic  are  presented  in  Chapter 
III.  The  notions  of  confounding  and  fractional  replication  which  are 
essential  in  the  logical  development  are  also  presented.  In  order  to 
present  factorial  experiments  in  which  the  factors  have  a  number  of 
levels  equal  to  the  power  of  a  prime  number  some  elementary  concepts 
of  Galois  field  theory  are  discussed. 

In  Chapter  IV  the  origin  and  structure  of  the  plans  given  in  the 
catalogue  are  presented.  The  efficiency  of  thi.  plans  is  described  and 
possibilities  of  blocking  are  discussed. 

The  construction  of  the  basic  plans  presented  in  the  catalogue  is 
described  in  Chapter  V.  Several  examples  of  orthogonal  main-eFect 
plans  constructed  from  the  basic  plans  are  given  and  an  index  of  J.e  plane 
which  can  be  obtai’'ed  from  the  catalogue  presented.  The  catalogue  of 
basic  orthogonal  main-effect  plans  then  conclude  the  report. 

D.  Notes  on  Terminology 

We  give  below  a  short  list  of  terms  which  occur  in  the  presentation 
with  some  explanation  of  their  meaning. 

(i)  A  Factor  designates  a  particular  force  which  is  varied  in  the 
total  investigation  at  the  will  and  under  the  control  of  the 
experimenter.  A  factor  is  also  called  an  input  variable  or  a 
controlled  variable. 
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(ii)  A  Quantitative  Factor  is  one  whose  values  can  be  arranged  in 

order  of  magnitude.  Such  values  can  usually  be  associated  with 
points  on  a  numerico.1  scale,  e.g.  temperatures  or  pressures. 
This  type  of  factor  is  also  called  a  continuous  factor  in  the 
literature. 

(iii)  A  Qualitative  Factor  is  one  whose  values  are  not  usually  arranged 
in  order  of  magnitude,  e.g.  type  of  dosage,  batches  of  material. 
Although  the  values  of  many  qualitative  factors  can  be  ordered 
according  to  a  particular  criterion  they  cannot  usually  be 
associated  with  points  cn  a  numerical  scale. 

(iv)  Levels  are  the  various  values  at  which  a  factor  is  examined,  e.g. 
the  levels  of  temperature  in  an  investigation  may  be  0°C,  50*^  C, 
100°  C  and  150°C. 

(v)  A  Treatment  Combination  is  one  of  the  possible  combinations  of 
levels  of  all  factors  under  investigation. 

(vi)  An  Experimental  Unit  is  that  entity  on  which  a  treatment  is 

applied.  In  experimentation  on  mice,  a  single  mouse  may  be  the 
unit.  In  agronomic  investigations  the  unit  is  frequently  a  plot  of 
land.  In  experimentation  on  a  chemical  process  the  unit  could  be 
the  system  for  a  prechosen  interval  of  time. 

{vii)  A  Trial  is  the  application  of  one  treatment  combination  on  one 
experimental  unit. 

(viii)  A  Response  is  the  result  a  trial  with  regard  to  a  particular 

attribute,  this  result  usually  being  expressed  numerically.  The 
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response  may  be  the  yield  of  a  process,  the  performance  of  a 
machine,  the  resistance  of  a  material  and  so  on.  Usually  there 
will  be  several  response  variables  for  each  trial. 

(ix)  An  Experiment  is  the  performance  of  a  planned  set  of  trials. 

(x)  A  Plan  is  a  set  of  treatnient  combinations. 

(xi)  The  Effects  of  a  factor  are  measures  of  the  change  in  response 
produced  by  a  change  in  the  level  of  the  factor.  When  a  factor  is 
examined  at  two  levels  only,  the  effect  is  the  difference  between 
the  average  response  of  all  trials  performed  at  the  first  level  of 
the  factor  and  that  of  all  trials  at  the  second  level.  If  there  are 
more  than  two  levels  the  differences  between  average  responses 
can  be  expressed  in  several  ways  e.  g.  linear  effects,  quadratic 
effects. 

(xii)  Error  is  the  variability  of  response  in  a  set  of  repetitions.  It 
usually  consists  of  components  of  different  origins,  e.  g.  failure 
of  units  to  be  identical,  failure  to  reproduce  treatment  combina¬ 
tions  exactly,  inaccuracies  of  measurement  of  responses. 
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II.  FITTING  LINEAR  MODELS  OR  REGRESSION  ANALYSIS 


The  basis  of  most  parametric  analyses  of  experiments  is  closely 


related  to  the  theory  of  fitting  linear  models  and  is  frequently  referred 
to  as  multiple  regression  .  Regression  analysis  can  be  defined  as  the 
estimation  or  predict!  jn  of  the  value  of  one  variable  from  the  values  of 
other  given  variables. 

The  assumption  m  regression  analysis  is  that  a  variable  y  may  be 


expressed  as  a  linear  function  of  some  known  variables  Xj,  x^. 


(which  may  be  functionally  related)  with  uncorrelated  random  deviations 

2 

which  are  distributed  around  zero  with  consvaiit  variance  cr  .  This 
linear  function  may  be  expressed  as 

y  =  3,  X,  +  a,x_  +  . . .  +  fl  X  +  e 
'a  *^1  1  ^2  2  '^p  p  a 

where  x, ,  x,,  . . . ,  x  take  on  a  particular  known  value  of  each  a  , 
say,  X  ,,  X  X  .  Frequently  x  ,  =  1  for  all  a. 

''  dl  (ip  Cll 

The  best  linear  unbiased  estimate  of  the  3’s  is  obtained  by 
minimizing  the  sum  of  squares  of  deviations 

.2 


S(y  -6,x,-S_x,-...-6x  ) 

al  ^2  a2  ^p  op 


The  term  regression  was  originally  introduced  to  describe,  partially, 
the  relationship  of  one  random  variable,  the  dependent  variable,  to 
another  random  variable,  the  independent  variable.  In  contexts  for  vhich 
regression  analysis  is  widely  used,  the  independent  variables  are  not 
random  variables,  so  the  term  is  not  entirely  appropriate. 
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This  procedure  is  known  as  che  method  of  least  squares.  Differentiating 


the  sum  of  squares  of  deviations  with  respect  to  each  of  the  p’s  in 
succession  the  follov/ing  equatinns  are  obtained; 


P,Sx  ,+p,Sx  ,x  +  p  Sx  ,x  =2y  x  . 

al  '^2  al  aZ  "^p  al  up  ai 

p,Sx  .X  ,  +  p,  Sx*",  +  ...  +pSx  ,x  =Sy  x 

al  aZ  ’^Z  aZ  aZ  ap  at. 


L  Sx  ,  X  +  p 
1  al  ap 


T  2  '■  o  X  + 

2  aZ  ap 


+  P  2x 


=  2y  X 
^a  ap 


These  equations  are  known  as  the  normal  equations.  If,  as  is  generally 

the  case  in  regression  problems,  the  x^^’s  aiv.  not  such  that  one  or  more 

linear  functions  of  them  are  zero,  then  a  unique  solution  of  the  above  set 

of  p  simultaneous  equations  exists.  In  order  to  solve  them,  first  solve 

p  sets  of  p  equations  the  first  set  of  which  is  written  as  follows,  using 

S..  =  S..  as  an  abbreviation  for  2x  .x  . 
ij  ji  ^  at  aj 


‘-1^11  ^2^12  ^ 


^1^12  ^2^22  +  •  •  •  +  %^2p  '  ° 


C,S,  +  C,S,  +  ...+CS_=0 

1  Ip  2  2p  p  pp 

Denote  the  solutions  of  these  equations  by  Cjj,  ...,  the 

first  subscript  indicating  that  this  is  the  solution  for  the  first  set  of 
equations  and  the  second  subscript  denoting  the  particular  C  solution. 
Next  solve  these  equations  with  unity  on  the  right-hand  side  of  the  second 
equation  and  zero  on  the  right-hand  side  of  all  the  other  equations,  the 
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solution  being  denoted  by  ^22’  '  *^2p*  Similarly  solve  the 

equations  with  unity  at  the  right-hand  side  of  the  third  eqiiation,  the 
fourth  equation,  and  so  or.  to  the  equation,  in  each  case  the  right- 
hand  side  of  all  other  equations  being  zero. 

The  solutions  can  be  arranged  in  a  pxp  square  as  follows: 


Cij  Cj2  •••  Cjp 
*"21  '^22***  ^2p 


C 

PP 


This  arrangement  of  the  solutions  is  known  as  a  matrix  and  is  the 

inverse  of  the  matrix  with  S. .  in  place  of  C.. .  The  solutions  for 

ij  ij 

P, ,  p,,  . . . ,  p  are  found  to  be  p.  =  S  C..  P.  where  P.  =  S  y  x  ... 

1  ^  P  i  j  ij  j  J  a  a  -.j 

It  will  be  noted  that  the  C..*s  are  derived  entirely  from  the  x..*s; 

that  is,  they  are  a  function  of  the  structure  of  the  observational  setup  and 

are  not  related  to  the  y’s  or  to  the  e*s.  The  quantities  estimating  the 

P’s  are  linear  functions  of  the  y  variables.  The  expectations  of  the 

A  A 

P’s  are  easily  found  to  be  the  corresponding  p’s,  the  variances  of  p.  to 

2  ^ 
be  and  the  covariance  of  any  two  p’s,  say  Pj^  and  p.  to  be 

?  Z 

C.j  <y  .  An  estimate  of  cr  is  derived  from  the  sum  of  squares  of 
deviations  about  the  estimated  values  and  is  given  by 

<r  =  s  =  — ^  S  (y  -  p,  X  ,  -  p,x  p  x  ) 

n-p  '’a  ol  aZ  '^p  up' 


=  5^  (Sy^  -  SpiPi) 
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where  2p.  is  the  sum  of  squares  removed  by  the  regression  on 
Xp  x^,  ....  Xp.  The  results  may  be  expressed  in  terms  of  the  analysis 
of  variance,  as  shown  in  Table  1. 

TABLE  1 

ANALYSIS  OF  VARIANCE 


Source 

m 

Sum  of  Squares 

Mean  Square 

Regxe  Jsion 

p 

P 

sp,p, 
i=l  *  ■ 

Zp.P./p  =  s^ 

Remainder 

(n-p) 

Difference 

2 

Difference /n-p  =  s 

Total 

n 

2 

2  y^ 
a=  1 

In  order  to  test  the  significance  of  the  regression  coefficients 

(the  p.  *s  )  the  random  deviations  e^  are  assumed  to  be  normally  and 

2 

independently  distributed  about  a  zero  mean  with  constant  variance  <r  . 

With  these  assumptions  the  significance  of  the  regression  coefficients 

can  be  tested  jointly  by  evaluating  the  mean  squares  in  the  analysis  of 

variance  and  comparing  the  ratio  s^/s^  to  the  F  distribution  with  p 

P 

and  (n-p)  degrees  of  freedom. 

With  the  extended  assumptions  on  the  random  deviations  e^  one  can 
also  construct  confidence  intervals  for  the  estimates  of  each  6. .  Since 
the  estimated  variance  of  is  C..  s  then  (p.  -  p)y  s  is 

distributed  as  Student's  t  distribution  with  (n-p)  degrees  of  freedom. 
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Hence  the  95%  confidence  intervals  on  are  given  by 


Pi  ±  t 


-  n-p,  95% 


SyT^.  . 
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Suppose  we  rename  the  regression  coefficients  p,,  p_,  . . . ,  p  , 

Pq+1>  •••.  Pp  and  we  %vish  to  test  whether  Pq+j.  Pq^.2.'  •  •  • .  Pp  could  be 
zero  making  no  assumptions  about  the  remaining  coefficients. 

The  procedure  is  as  follows: 


(i)  Estimate  the  regression  coefficients  in  the  model 


pjXj  +  p,x 


2''2 


+  P  X  +  e 
P 


obtaining  ^2‘ 
the  regression  on 


The  sum  of  squares  removed  by  the 

P  '• 

. . ,  X  IS  equal  to  S  P.  P.  . 

P  i=l  ^  ^ 


(ii)  Estimate  the  regression  coefficients  in  the  model 


~  ^1  ^  P2^2  '*'•••  Pq^q  ^ 


obtaining  p^,  p| ,  , . ,  ,  p*  . 
regression  on  Xj,  x^, 


The  sum  of  squares  removed  by  the 

X  is  equal  to  B  0*P  . 

q  ..11 

^  1=  1 


(iii)  Construct  the  analysis  of  variance  given  in  Table  2. 
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TABLE  2 


ANALYSIS  OF  VARIANCE 


Source 

d.  f. 

S\im  of  Squares 

Mean  Square 

Regression  on  x,,...,x^ 

q 

2 

q 

8 

q 

Regression  on  "'^p 

p  q 

2 

after  fitting  Xj, . . . ,  x^ 

P-q 

1=1  1=1 

®d 

p 

2 

Regression  on  Xj,...,x^ 

P 

s  p.p. 

•  it  t 

1=1 

s 

p 

Remainder 

n-p 

Difference 

2 

s 

Total 

n 

2 

S  y 

1  P 
o- 1 

To  test  the  hypothesis  that  •  •  • »  P  zero  we  utilize  the 

2  2 

fact  that  under  the  hypothesis  that  they  are  zero  the  ratio  s^/s  will  be 

distributed  as  F  with  (p-q)  and  (n-p)  degrees  of  freedom  and  thus 
2  Z 

compare  s^/s  with  the  v<iiue  in  the  F  table  corresponding  to  (p-q) 
and  (n-p)  degrees  of  freedom. 

the  usual  regression  test  devised  to  test  whether  deviations  about 
the  mean  have  a  regression  on  the  independent  x  variates  may  be 
deduced  from  the  above  discussion.  The  complete  hypothesis  is  that 


0  X 

P 


+  e 

a 


22 


and  the  restricted  hypothesis  that 


y  =  8,  X  ,  +  e 
'^1  al 

where  is  unity  for  all  values  of  a-  The  estimate  is  y  and  the 
sum  of  squares  due  to  the  regression  on  x^  (i.  e.  the  sum  of  squares 
due  to  the  mean)  is  y  ZJy .  The  "correction  for  the  mean"  y  Sy ,  with 
one  degree  of  freedom  may  be  deducted  from  the  total  sum  of  squares 
and  the  analysis  is  given  in  Table  3. 

TABLE  3 

ANALYSIS  OF  VARIANCE 


Source 

B 

Mean  Square 

p  * 

? 

Regression  on  x*, .  • ,  ,x 
c  p 

P-1 

2  Pi  Pic 
i^2  ^ 

s 

r 

Remainder 

n-p 

Difference 

2 

s 

Total 

n-1 

n  2  _  n 

S  y  -  y  2  y 
a= 1  a= 1 

P.  =  P.  -  X.  Sy  denotes  the  sum  of  products  around  the  mean, 
ic  i  11  a 

Much  of  the  proceeding  discussion  can  be  simplified  through  the  use 
of  matrix  notation.  The  linear  function  expressing  y  as  a  function  of 
the  X  variates  may  be  written  as 

y  =  Xp  +  e 
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where 


'yi' 

^11 

Xi2  -  . . 

X, 

Ip 

'Pi' 

yz 

• 

-  .  . 

•  • 

P2 

y  = 

ya 

,  x  = 

\2  •  •  • 

X 

ap 

,  . 

.  P  = 

• 

^nl 

^n2  "  • 

L _ 

1 

o, 

-  CO. 

_ 

The  sum  ot  square?  to  be  minimized  is 


e»e  =  (y  -  xp)»  (y  -  Xp)  =  y*y  -  2p*X‘y  +  p*X»Xp. 


The  normal  equations  are 

.Sp  =  X»y 

where  S  =  X*X. 

If  S  is  non-singular  then  p  =  S'’^X*y  and  the  variance-covariance 

2-1  ^ 

matrix  of  the  estimates  is  equal  to  a-  S  .  The  estimate  of  o-  is 
given  by 

0-“  =  (y  -  Xp)*(y  -  Xp)/(n-p)  =  (y’y  -  P’ X*y}/{>i-P)  • 

The  results  presented  above  will  be  utilized  in  the  chapters  which 
follow  and  matrix  notation  will  be  used  whenever  it  simplifies  the 
presentation. 

2 

There  arises  the  problem  of  how  should  c  be  estimated  if  n  =  p 

(i.  e.  the  number  of  parameters  to  be  estimated  equals  the  number  of 

2 

trials).  The  estimate  of  v  based  upon  the  s\im  of  squares  of 

deviations  about  the  estimated  values  of  the  parameters  is  sometimes 

■*■2 

called  pure  error.  If  n  =  p  it  is  clear  from  the  formula  for  <r  that  no 
estimate  of  pure  error  can  be  derived  from  the  experiment..  In  such  a 
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situation  there  ar^;  two  possible  ways  of  resolving  the  problem.  First, 
the  experimenter  may  be  Investigating  a  process  for  which  the  experi¬ 
mental  error  is  already  knovm.  In  this  case  the  error  obtained  from 
prior  information  may  be  used  as  an  estimate  of  a-  .  This  estimate  of 
QT^  can  then  be  based  on  infinite  degrees  of  freedom  and  the  estimation 
and  test  of  significance  procedure  can  be  made  as  if  the  estimate  of 
experimental  error  had  been  obtained  from  the  experiment  itself. 

A  graphical  procedure  fo.  analysing  factorial  experiments  developed 
by  Daniel  (1959)  may  be  u.seful  in  obtaining  a  rough  estimate  of  error. 

This  procedure  uses  a  half-normal  grid  on  which  to  plot  the  absolute 
values  of  the  contra.'ts  defining  the  ma.-a  effects  and  interactions  of  a 
factorial  experiment.  If  these  contrasts  are  arranged  in  order  of 
absolute  magnitude  and  plotted  on  half-normal  probability  paper  they  should 
fall  along  a  straight  line,  if  all  factoi-s  have  no  effects. 

A  half-normal  grid  can  be  prepared  by  taking  a  sheet  of  arithmetic 
(normal)  probability  paper,  deleting  the  printed  probability  scale  P,  for 
the  range  P  <  50%  and  replacing  each  value  of  P  >  50%  by  the  corre¬ 
sponding  value  of  P*  =  2P  -  100.  The  relatio.n 

P*  =  (i  -  |)/N:  i  =  i,  2,  . .  .,  N, 

where  N  is  the  number  of  main  effects  and  intera.rtions  to  be  estimated, 
is  used  for  plotting  the  empirical  distribution  oi  contrasts.  The  abscissae 
are  the  absolute  values  of  the  contrasts. 

Under  the  null  hypothesis  that  all  factors  have  no  effects  the  standard 
error  of  each  contrast  ,  could  be  roughly  estimated  by  Ihe  contrast  for 
which  P’  is  most  nearly  0.  683.  If  it  is  known  that  some  effects  or 
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interactions  are  likely  to  be  real  and  it  appears  from  the  graph  that  they 
are,  then  they  should  be  judged  real  and  the  remaining  contrasts  used  to 
determine  the  standard  error  i.  e.  reduce  N  by  the  number  of  real 
effects  and/or  interactions  in  the  formula  P’  =  (i  -  ■^)/N  and  then 
estimate  the  standard  error  of  a  contrast  by  the  absolute  value  of  the 
contrast  for  which  the  new  P*  is  most  nearly  0.  683.  If  a  straight  line 
is  drawn  through  the  origin  and  the  absolute  value  of  the  contrast  for 
which  P*  is  most  nearly  o.  683  a  rough  idea  of  which  effects  and 
interactions  are  significantly  large  can  be  obtained.  These  will  fall  far 
to  the  right  of  the  line.  An  estimate  of  the  experimental  error,  (r^,  can 
be  obtained  from  the  formula 

*'7. 

=  (r^/{N+  1)  . 

The  estimate  of  experimental  error  is  based  upon  N  degrees  of  f.eedom 
and  although  it  is  approximate  and  deduced  by  subjective  reasoning,  it 
does  give  some  inf<  rmation  about  the  experiment  that  would  not  be 
forthcoming  without  an  estimate  of  error.  The  reader  who  is  interested  in 
this  technique  can  find  many  illustrative  examples  of  its  use  in  the  paper 
by  Daniel  (1959). 
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III.  FACTORIAL  EXPERIMENTS 


When  an  experiment  involves  several  factors,  the  effects  of  all 
factors  on  a  characteristic  of  interest  may  be  investigated  simultaneously 
by  varying  each  factor  so  that  all  or  a  suitable  subset  of  all  possible 
'■ombinations  of  the  factors  are  considered.  An  experiment  in  which  this 
procedure  is  used  is  known  as  a  factorial  experiment. 

A.  Factorial  Experiments  with  Factors  at  Two  Levels 

The  simplest  and  most  common  factorial  experiments  involve  factors 
which  occur  at  two  levels.  The  two  levels  of  a  factor,  may  be  denoted  by 
0  and  1.  A  treatment  is  denoted  by  a  particular  combination  of  levels, 
one  level  from  each  factor.  The  treatment  combination  for  which  all  the 
factors  occur  at  the  0  level  can  be  simply  denoted  by  -{1).  The  1  level 
of  a  factor,  say  factor  A,  can  also  be  represented  by  the  lower  case 
letter  a.  A  factorial  experiment  involving  three  factors  A,  B  and  C 
each  at  two  levels  would  consist  of  the  following  treatment  combinations: 
(1),  a,  b,  ab,  c,  ac,  be  and  abc.  In  these  combinations  the  presence  of 
a  letter  indicates  that  the  corresponding  factor  occurs  at  its  1  level  and 
the  absence  of  a  letter  indicates  the  corresponding  factor  occurs  at  its  0 
level. 

The  main  effect  of  factor  A  is  defined  to  be  the  t’ifference  between 
the  mean  of  the  yields  at  the  1  level  of  factor  A  and  the  mean  of  the 
yields  at  the  0  level  of  factor  A. 

Hence  the  main  effect  of  A  is  ^(a  +  abt  ac  +  abc)  -  i(  (1)  +  b+ c  +  be ) 
which  can  also  be  written  as 
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A  =  i{a-l)(b+l){c+l) 

where  the  expression  is  to  be  expanded  algebraically  and  the  responses 

substituted  for  the  treatment  symbols.  The  effects  and  interactions  of 
3 

the  2  factorial  experiment  are  given  by 

A  =  j(a-l)(b+l){c+l) 

B  =  i(a+l)(b-l)(c+l) 

AB  =  i(a-l)(b-lKc+l) 

C  =  i{a+l){b-rl)(c-l) 

AC  =  i(a-l){b+l)(c-l) 

BC  =  i(a+l)(b-l){c-l) 

ABC  =  i(a-l)(b-l)(c-l) 

a  minus  sign  appearing  in  any  factor  on  the  right  if  the  letter  is  present 
on  the  left.  We  will  adhere  to  the  convention  that  treatment  combinations 
are  represented  by  lower-case  letters  and  effects  and  interactions  by 
capitals. 

It  will  be  noted  that  the  effects  and  interactions  are  seven  mutually 
orthogonal  contrasts  of  the  responses  of  the  eight  treatment  combinations. 


28 


/ 


e 


(X)  a  b  ab  c  ac  be  abc 

4A  -  +  -  +  -  + 

4B  --  +  +  --  +  + 

4AB  +  _-  +  +  --  i- 

4C  -_--  +  +  +  + 

4  AC  +-+--+-+ 

4BC  +  +  --  -_  +  + 

4ABC  -■  +  -  +  --■!■ 

Orthogonality  of  two  linear  contrasts  may  be  defined  as  follows: 

Consider  two  linear  functions,  C^  and  of  the  variates 

X, ,  Xt,  . . . ,  X  where  the  x*s  have  the  same  variance  and  are 
la  n 

uncorrelated. 

C,  =  a,  X,  +  tto  X-  +  . . .  +  a  X 

X  X  1  b  b  11  n 

^2  ~  Pi  ^1  ^2  ^n  ^n 

where  and  may  assume  any  values,  not  all  zero.  A  necessary 
and  sufficient  condition  that  the  two  linear  fur-.ctions  be  orthogonal  is 

n 

S  o.  p.  =  0 
i=l  ^  ‘ 

If  the  mean  response  of  the  eight  treatment  combinations  is  denoted 
by  p.  the  effects  and  interactions  are  represented  by 


8  p 

1 

1 

1 

1 

1 

1 

1 

1  " 

r«i 

4  A 

-1 

1 

-1 

1 

-1 

1 

-1 

1 

a 

4  3 

-1 

-1 

1 

1 

-1 

-1 

1 

1 

b 

4  AE 

- 

1 

—  1 

I 

—  i 

1 

1 

-1 

-1 

1 

ab 

4  C 

-  i 

-1 

-1 

t 

“  a 

1 

1 

1 

1 

c 

4  AC 

1 

-1 

1 

-1 

-1 

1 

-1 

1 

ac 

4  DC 

1 

1 

-1 

-1 

-1 

-1 

] 

1 

be 

_4ABC 

-1 

1 

1 

-1 

1 

-1 

-1 

1 

abc 

V/ith  ii  factors  A,  B,  C,  JJ,  etc.  the  effects  and  interactions  may  be 
represented  by 

X  =  (a  +  l)(b  +  l)(c  +  T)(d  +  1)  . . . 

-  -  -  - 

where  the  sign  in  each  bracket  is  positive  if  the  corresponding  cap'-.a! 
letter  is  not  contained  in  X  and  negative  if  it  is  contained  in  X,  and  the 
whole  expression  or  the  right-hand  side  is  to  be  expanded  algebraically 
and  the  yields  substituted  in  place  of  the  corresponding  treatment 
combinations. 

The  choice  of  the  divisor  in  the  above  expression  is  a  matter  of 
convention  only  and  depends  upon  the  definition  of  an  effect  or  interaction. 
Here  we  have  defined  an  effect  or  interaction  on  the  basis  of  the  difference 
letween  two  experimeilVal  units. 

The  response  of  a  treatment  combination  may  be  written  as 
a;  bj  Cj,  . . .  where  absence  is  denoted  by  the  subscript  taking  the  value 
zero  and  presence  by  the  subscript  taking  the  value  unity.  Then 
a.b^Cj^  ...  =  p  +  |a  +  t^^B  +  2^  il-AC  +  ... 
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where  the  sign  on 

4a 

IS 

-  if 

i  =  0 

and 

if 

i  =  1 

on 

is 

-  if 

j  =  0 

and  t 

if 

j  =  1 

on 

IS 

-  if 

k  =  0 

and  -1 

if 

k  =  1 

and 

SO 

on^ 

and  the  sign  on  a  term  involving  several  letters  is  the  product  of  the 
signs  on  the  individual  letters. 

If  a  2”  experiment  is  replicated  r  times  in  randomized  blocks  of 

2*'  plots  each  effect  or  interaction  is  estimated  by  the  mean  of  r  2  ~ 

11“  1 

responses  minus  the  mean  of  r  2  responstr  and  therefore  has  a 

-2 


variance  of 


f  I  .  1  \  2  2  ,  .,n- 

/ - }  +  - ;  jo-  =0-  /r  2 

^2“-^  r2”“7 


Furthermore  the 


estimates  of  effects  and  interactions  are  uncorrelated  so  that  the 
variance  of  any  linear  function  of  them  can  be  easily  obtained. 


B.  Factorial  Experiments  with  Factors  at  More  than  Two  Levels 
The  3^  system 

With  factors  at  three  levels  the  effect  of  any  one  factor  may  be 
expressed  in  several  ways.  First  the  response  at  each  level,  where  the 
level  is  represented  by  0,  1  or  2,  can  be  expressed  as  a  deviation  from 
the  mean  response  at  the  three  levels,  giving  say  Aq,  Aj  and  where 
Aq  Aj  -t  =  0.,  Another  approach  is  that  the  main  effect  of  a  factor 
can  be  represented  by  independent  comparisons  among  the  means  corre¬ 
sponding  to  the  different  levels  of  the  factors.  Among  three  independent 
quantities  there  are  two  independent  comparisons.  The  comparisons 
which  are  of  interest  will  depend  upon  the  nature  of  the  factors,  in 
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particular  whether  they  are  qualitative  or  quantitative. 

If  the  levels  are  qualitative  and  the  0  level  denotes  the  control  and 
the  other  two  denote  treatments  the  comparisons  of  interest  may  be 
(i)  a  comparison  of  the  two  treatments  and  (ii)  a  comparison  of  the 
average  of  the  two  treatments  with  the  control.  These  comparisons  can 
be  expressed  as 

A*  =  2(2ao-a^  -  a^) 

A«»  =  aj  -  a, 

respectively,  where  a^,  a^  and  a^  denote  the  responses  with  factor 
A  at  the  0,  1  and  7.  levels. 

For  most  quantitative  factors  the  ccmparisons  of  interest  will  be 
those  giving  the  most  information  on  the  relation  between  the  responses 
and  the  levels,  namely  the  slope  and  the  curvature.  This  can  be 
represented  by  a  polynomial  expression  y  =  +  a  j  x  +  ,  where  x 

denotes  the  levels  of  the  factor  and  y  is  the  response  variable.  The 
linear  and  quadratic  effects  of  factor  A  may  be  written  as: 

■^Q  “  ^^2  "  ^®1  ^0 ^  ■ 

The  quadratic  effect  is  the  linear  contrast  among  a^,  a^  and  which 
is  orthogonal  to  the  linear  effect. 

Now  consider  two  quantitative  factors  A  and  B,  each  at  three 
equally  spaced  levels.  The  interaction  of  these  two  factors  will  be  the 
interaction  of  a  3x3  table  and  will  have  four  degrees  of  freedom. 


These  four  degrees  of  freedom  may  be  separated  into  orthogonal 
contrasts  each  with  a  single  degree  of  freedom. 

■  V 

=  (^2  -  ^0)0^2  “  +  bjj) 

=  (a2-2aj  +  ay)(b2-2b^+bQ) 

This  system  of  expressing  .ne  results  may  be  extended  indefinitely. 

Several  conventions  have  been  used  to  define  the  main  effects  and 
interactions,  each  convention  having  some  merit.  One  common 
convention  adopted  is  to  define  the  effects  and  interactions  on  the  basis  of 
the  difference  between  two  experimental  units.  Adopting  this  convention 
the  main  effects  and  interactions  of  an  experiment  on  two  three-level 
lactors  A  and  B,  are  given  by 

Aj^  =  ^  (a^  “  ®Q)(bQ  +  bj  +  b^) 

B^,  —  ^  (s-Q  ^  "  ^0^ 

®Q  "  6  ^  ^1  ^  ^z'^^O  "  ^^1  ^ 

■^L^L  ”  2  ^^2  ”  ®0^^^2  ”  ^0^ 

^Q®L  =  +  V 

The  convention  adopted  does  not  alter  any  tests  of  significance 
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performed  on  the  parameters  and  therefore  need  be  of  little  concern  to 
the  experimenter. 

There  is  a  class  of  experiments  involving  quantitative  and  qualitative 
factors  in  which  the  treatment  combinations  have  an  appearance  of 
consisting  of  a  full  set  of  factorial  combinations  but  are  not  in  fact  so. 

A  simple  example  Ox  3  type  is  that  in  which  there  are  three  eoually 
spaced  amounts,  including  a  zero  amount  of  a  particular  treatment 
administered  by  three  me'-iods.  Since  the  zero  amounts  of  the  treatment 
administered  by  the  three  methods  are  identical  treatments  there  are  only 
seven  different  treatment  combinations  and  not  nine.  The  experimenter 
must  consider  whether  he  .should  use  the  nine  treatment  combinations  as 
though  they  were  all  distinct  or  only  the  set-cn  distinct  combinations,  and 
further  he  should  consider  the  method  of  analysis  in  each  case.  For  a 
more  detailed  discussion  of  this  type  of  experiment  the  reader  is 
referred  to  section  18.8  of  Kempthorne  (1952), 

We  now  present  a  formal  method  of  defining  effects  and  interactions. 
Consider  the  case  of  three  factors  A,  B  and  C  each  at  two  levels  0 
and  1.  The  eight  treatment  combinations  (1),  a,  b,  ab,  c,  ac,  be,  abc 
may  be  represented  by  the  points  (0,0,0),  (1,0,0),  (0,1,0),  (1,1,0), 
(0,0,1),  (1,0,1,,  (0,1,1)  and  (1,1,1)  respectively,  in  Euclidean 
space  with  axes  Xj,  x^  and  x^,  the  first  coordinate  referring  to  the 
level  of  factor  A,  the  second  to  the  level  of  factor  B  and  the  third  to  the 
level  of  factor  C.  The  effects  and  interactions  defined  previously  have  a 
simple  algebraic  interpretation.  The  effect  of  A  is  the  comparison  of  the 
treatment  combinations  for  which  Xj  =  0  with  those  for  which  Xj  =  1. 
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Likewise  the  effect  of  B  is  the  comparison  of  the  treatment  combinations 
for  which  =  0  with  those  for  which  x_  =  1  and  the  effect  of  C  is  the 
comparison  ol  the  treatment  combinations  for  which  x.  =  0  with  those 
for  which  -  i.  The  interaction  AB,  for  example,  is  in  the  former 
notation  the  comparison  among  treatment  combinations, 

(1)  +  c  +  j.b  +  abc  -  a  -  h  -  ac  -  be 

i.  e.  of  the  points  (0,0,0),  (0,0,1),  (1,1,0)  and  (1,  1,  1)  versus  the 
points  (LO.O),  (0,1,0),  (1,0,1)  and  (0,1,  l)i.  For  the  points  (0,0,0) 
and  (0,0,1),  x^  +  X2  =  0  and  for  the  points  (1,1,0)  and  (1,1,1), 

X,  +  x^  =  2  and  for  the  other  four  points,  =  1.  If  the  numbers 

are  reduced  modulo  2,  that  is,  any  number  is  replaced  by  the  remainder 
when  it  is  divided  by  2,  the  interaction  is  the  comparison  of  those 
treatment  combinations  for  which  Xj  +  x^  =  0  (mod  2)  versus  thos-  for 
which  Xj  +  x^  =  1  (mod  2).  It  is  easily  verified  that  the  effects  and 
interactions  are  bas  d  on  a  comparison  of  two  groups  of  treatment 
combinations  given  by  the  equations  in  Table  4. 
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TABLE  4 


EQUATIONS  REPRESENTING  EFFECTS  AND  INTERACTIONS 


Effect  or 
Interaction 

Left-Hand  Side 
of  Equation 

A 

^1 

B 

^2 

AE 

Xi  +  X2 

C 

^3 

AC 

Xi  +  x^ 

BC 

^2  ^3 

ABC 

Xj  +  X2  -i-  Xj 

For  example  the  treatment  combinations  entering  ABC  with  a  minus 
sign  are  (1),  ab,  ac  and  be  and  for  these  x^  +  X2  +  =  0  (mod  2)  and 

the  treatment  combinations  entering  with  a  plus  sign  are  a,  b,  c,  and 
abc  for  which  +  x^  +  x^  =  1  (mod  2). 

The  above  approach  for  the  system  suggests  the  appropriate 
approach  for  the  3*^  system.  Consider  the  arrangement  of  the  nine 
treatment  combinations  with  two  factors  each  at  three  levels. 


Level  of  factor  A 

(0,0) 

(1.0) 

(2,0) 

Level  of 
factor  B 

(0,1) 

(1.1) 

(2,  1) 

(0.2) 

(1.2) 

(2.  2) 

4 

^2 
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The  main  effect  of  factor  A  can  be  represented  by  the  comparisons 
among  three  means;  those  for  which  =  0,  for  which  =  1  and  for 
which  =  2.  A  representation  of  these  effects  as  two  linearly 
independent  numbers  may  be  obtained  by  considering  each  mean  as  a 
deviation  from  the  over-all  mean.  The  interaction  of  factors  A  and  B 
has  four  degrees  of  freedom.  These  four  degrees  of  freedom  can  be 
considered  from  the  point  of  view  of  the  completely  orthogonalized  5xi 
square; 


The  comparisons  among  the  columns  give  the  effect  of  factor  A,  and 

among  the  rows  the  effect  of  factor  B.  Those  among  the  Latin  letters 

and  those  among  the  Greek  letters  each  with  two  degrees  of  freedom 

represent  the  four  degrees  of  freedom  for  the  interaction  of  the  two 

factors.  Consider  the  following  grouping  given  by  the  Latin  letters;, 

(0,0L  (2,1),  (1,2)^  versus  (1,0),  (0,1^,  (2,2)  versus  (2,0),  (0,2),  (l.i). 

For  this  grouping  the  comparisons  are  among  those  treatment  combinations 

for  which  Xj  +  x^  =  0,  =  1,  =  2  (mod  3).  Similarly  the  comparisons  among 

the  Greek  letters  are  comparisons  among  the  treatment  combinations  for 

which  Xj  T  2X2  =  0,  =  1,  =  2  (mod  3). 

The  pair  of  degrees  of  freedom  corresponding  to  the  equations 

Xj  +  x^  =  0,  =  1,  =  2  may  be  denoted  by  the  symbol  AB  and  the  pair 

2 

corresponding  to  Xj<-2x2=  0,  =1,  =2  by  AB  .  The  interaction  degrees 
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of  freedom  may  also  be  represented  by  BA  and  BA  respectively.. 

It  is  easily  verified  that  the  comparisons  among  the  groups  of  treatment 

2 

combinations  represented  by  BA  and  BA  are  the  same  as  those 

2 

represented  by  AB  and  AB  respectively.  It  is  necessary,  in  order 
to  obtain  a  complete  and  unique  enumeration  of  the  pairs  of  degrees  of 
freedom,  to  adopt  the  rule  that  an  order  of  the  letters  is  to  be  chosen  in 
advance  and  i.hat  the  power  of  the  first  letter  in  a  symbol  must  be  unity. 

K  the  power  of  the  first  le'.er  of  a  symbol  is  2  then  by  squaring  the 
symbol  and  using  the  rule  that  any  letter  cubed  is  to  be  replaced  by  unity 
the  power  of  the  first  letter  will  be  unity.  This  process  may  be  extended 
indefinitely.  For  three  factors  the  results  are  shown  in  Table  5. 

The  extensions  are  quite  straightforward  and  need  not  be  enumerated. 
For  the  system  there  are  n  independent  factors  and  their  generalized 
interactions,  giving  rise  to  (3*^  -l)/2  symbols  each  representing  two 
degrees  of  freedom.. 

The  symbols  used  above  to  denote  pairs  of  degrees  of  freedom  can 
also  be  used  to  denote  the  magnitudes  of  effects  and  interactions.  Each 
symbol  represents  a  comparison  among  three  groups  of  3^”^  treatment 
combinations,  examples  of  which  are: 
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TABLE  5 


EQUATIONS  REPRESENTING  EFFECTS  AND  INTERACTIONS 


Effect  or 
Interaction 

Left-Hand  Side 
of  Equation 

A 

^1 

B 

^2 

AB 

X1  +  X2 

AB^ 

Xi  +  2x2 

C 

^3 

AC 

Xj  +X3 

AC^ 

Xj  +  Zx^ 

BC 

X2  +  X3 

BC^ 

x^  +  Zxj 

ABC 

Xi  +  X2  +  X3 

ABC^ 

Xj  +  x^  +  2x, 

AB^C 

Xj  +  Zx^  +  X3 

AB^C^ 

Xj  +  Zx^  +  2X3 

Aq  =  (mean  of  treatment  combinations  for  which  Xj  0  (mod  3)  ) 
-  (mean  of  all  treatment  combinations) 

ABq  =  (mean  of  treatment  combinations  for  which  x^  +  x^  =  0 
(mod  3)  )  -  (mean  of  all  treatment  combinations) 

2 

ABj  =  (mean  of  treatment  combinations  for  which  Xj  +  Zx^  =  1 
(mod  3)  )  -  (mean  of  all  treatment  combinations) 
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AB^C-,  =  (mean  of  treatment  combinations  for  which  +  +  “  2 

(mod  3)  )  -  (mean  of  all  treatment  combinations). 

With  these  definitions  the  response  of  treatment  combination 

a.  b.  c,  in  terms  of  effects  and  interactions  is 
1  j  k 

a.b.c,  =  U  +  A.  +  B.+  AB. ,  .  +  AB^t  +  C,  +  AC.,,  +  AC^,.  +  BC . , , 

i  j  k  1  j  i+j  i.+  2j  k  i+k  i+2k  j+k 


+  BC?^,,  +  ABC.^.^,  +  ABC?,.^,,  +  AB^C..,.., 
jtZk  i+1+k  :-^j-i~2k  i+Zj+k 


where  all  subscripts  are  reduced  modulo  3  and  p  is  the  mean  of  all 
combinations.  For  example,  tne  response  of  treatment  combination 
^1  bfl  c^  is  given  by 

ai  bo  c^  =  p  +  A  J  +  Bq  +  ABj  +  AB  ^  +  C2  +  ACq  -t  AC^  +  BC2  + 

+  ABCq  +  ABC2  +  AB^Cq  +  AB^C^  . 

Thus  it  is  possible  to  express  any  linear  contrast  of  the  responses  in 

terms  of  the  effects  and  interactions. 

Now  suppose  that  the  treatment  combinations  are  tested  the  same 

number  of  times  in  a  randomized  block  trial.  Then,  with  an  additive 

model,  the  observed  response  will  be  equal  to  a  true  response  plus  an 

error.  T’ne  errors  may  be  regarded  as  uncorrelated  with  mean  zero  and 
2 

constant  variance  cr  .  Then  the  best  estimate  of  any  contrast  of  the 
true  responses  is  the  same  contrast  of  the  observed  means. 

The  only  estimable  functions  of  the  parameters  are  functions  of  the 
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type  where  a  is  one  of  the  set  of  symbols  A,  B,  AB,  AB^, 


ABC  etc.  and  the  i,j  have  values  equal  to  0,  1  or  Z.  It  is  easily 
verified  that  the  estimates  of  quantities  a.  -  and  -  P^,  where 
a,p  are  different  ones  of  the  set  of  symbols  are  uncorrelated. 


2 

Consider  the  nine  treatment  combinations  of  the  3  factorial 


experiment  written  in  terms  of  effects  and  interactions. 

2 

Uq  bg  =  p.  +  Aq  +  +  ABq  +  .A.Bq 

2 

aQ  bj  -  p  +  Aq  +  Bj  -S-  AB^  +  AB^ 

2 

Uq  b^  =  p  +  Aq  +  B2  +  AB^  ABj 

aj  bQ  =  p  +  Aj  +  Bq  +  ABj  AB^ 

Ui  bj  =  p  +  Aj  +  Bj  +  AB2  +  ABq 

2 

a.  :=  +  Aj  +  +  AB^  +  AB^ 

2 

b^  -  \i  +  A^  +  Bq  +  AB2  +  -^^2 

82  bf  =  p  +  A2  +  Bj  +  ABq  +  ABj 

2 

^2  ^2  ~  ^  "^2  ^  ^2  ^  ABj  +  ABq 

The  estimate  of  A2  ■  Aq,  say,  is  clearly  equal  to  a  constant  limes 

(^2  ^0  ^2  ^2  ^2  "  ^0  *^0  ■  ^0  ■  ^0  ^2^ 
and  the  estimate  of  B2  ”  equal  to  a  constant  times 

(ao  b2  +  aj  b2  +  a2  b2  -  aQ  b^  -  a^  b^  -  82  b^) 

The  coefficients  of  the  treatments  for  the  above  two  contrasts  are, 
apart  from  the  constant  multiplier 

4t 


^0*^0 

^0*^1 

^0^2 

^1*^0 

ajb. 

^1^2 

a^bo 

^2*^2 

■^2"‘^0 

-1 

-1 

0 

0 

0 

1 

1 

1 

B,-Bj 

0 

-1 

t 

0 

-  i 

1 

0 

-1 

1 

Since  the  sum  of  the  products  of  corresponding  coefficients  is  zero 
the  two  contrasts  are  orthogonal. 

Among  the  three  deviations  from  the  overall  mean  Aq,  A^  and  A^, 
say,  there  are  two  independent  contrasts.  These  can  be  represented  by 
the  contrasts  A^  -  Aq  and  A^  -  2Aj  +  Aq.  If  the  same  types  of 
contrasts  are  utilized  for  each  of  the  other  symbols,  it  is  possible  to 
obtain  eight  orthogonal  contrasts. 

Since  in  each  of  the  two  systems  of  defining  effects  and  interactions 

the  treatment  combinations  can  be  written  in  terms  of  the  effects  and 

interactions  it  is  not  difficult  to  determine  the  relationship  of  one  /stem 

2 

to  another.  For  example,  in  the  3  experiment  on  factors  A  and  B 

A.  =  a,  b„  +  a,  b-  +  a-  b*  -  a-  b,,  -  a,,  b,  -  a„  b- 

L  rO  21  22  00  01  02 

and  A^-AqS  a^  +  a^  bj  +  a^  b^  -  a^  b^  -  a^  b^  -  a^  b^ 

Thuo  Aj^  =  A^  -  Aq.  Similarly,  it  can  be  shown  that  Aq  =  Aq  -  2Aj 
Now,  =  a^  b^  +  b^  -  a^  b^  -  a^  b^.  If  these  four  treatment 

combinations  are  substituted  in  the  equation 

a.  b.  =  11  +  A.  +  B.  +  AB.  .  +  AB. , 

1  J  1  J  i+J  i+2j 

it  is  easily  demonstrated  that 
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The  following  presentation  is  a  straightforward  generalization  of  the 
2*^  and  3^  systems.  The  generalization  from  the  3*^  system  to  the 
system,  where  p  is  a  prime  number,  can  be  seen  fairly  easily, 
without  introducing  proofs.  The  proofs  are  based  on  the  properties  of 
Galois  fields  which  will  be  given  later. 

Represent  the  treatment  combination  by  numbers  x,  x-  .  .  .  x  , 

1m  h 

th 

where  is  the  level  oi  the  i  factor  in  the  particular  combination. 
The  numbers  take  on  values  from  0  to  (p-1).  All  the  numbers  arc 
reduced  modulo  p,  that  is,  a  number  greater  than  (p-l)  is  replaced  by 
the  remainder  after  division  by  p.  The  (p^-1)  degrees  of  freedom 
among  the  p^  treatment  combinations  may  be  partitioned  into 
{p^-l)/(p-l)  sets  of  (p-l)  degrees  of  freedom.  Each  set  of  (p-l) 
degrees  of  freedom  is  given  by  the  contrasts  among  the  p  sets  of  , 


p  equations: 


(mod  p) 


The  Uj,  ’s  must  be  positive  integers  between  0  and  (p-l),  not  all 
c’cal  to  zero  and  for  uniqueness  the  coefficient  of  the  first  a.  that  is 
not  zero  equals  unity. 

Two  sets  of  (p-l)  degrees  of  freedom  resulting  from  equations  with 
left-hand  sides  Sa.  ::  and  will  be  orthogonal  unless  p.  -  ka. 

for  each  i  .  This  can  easily  be  seen  because  the  two  equations, 
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z^a.  X.  =  k 

^  (mod  p) 

S6.  X  =  m 

1 

n  2 

will  be  satisfied  by  p  treatment  combinations,  if  p.  is  not  equal  to 
a  constant  multiplier  of  a.  . 

*^1  °’2  ®n 

The  symbol  A  B  . .  .  K  ,  which  corresponds  to  the  equations 

whose  left-hand  side  is 

a,  X,  +  a,  X,  +  . .  .  +  a  x  , 

112  2  n  n 

denotes  a  set  of  (p-1)  degrees  of  freedom,  the  power  of  the  first  letter 
occurring  being  restricted  to  be  unity. 


Galois  field  theory 

In  order  to  obtain  a  procedure  for  investigating  factors,  each  having 
s  levels,  where  s  -  ,  a  knowledge  of  group  theory  is  essential. 

A  set  of  s  elements  vt^,  u^,  . .  .  ,  ^  is  said  to  be  a  finite  field  of 

order  s  if  the  following  properties  hold: 

(i)  The  set  is  closed  under  addition  and  multiplication,  i.  e.  if  u. 
and  u.  belong  to  the  set  then  so  do  u.  +  u.  and  u.  u. . 

(ii)  Addition  and  multiplication  are  commutative,  i.  e. 

u.  +  u.  -  u.  +  u.  and  u.  u.  =  u.  u.  . 

t  J  J  I  1  j  1  I 

(iii)  Addition  and  multiplication  are  associative,  i.  e. 


(iv)  The  distributiv"  law  holds,  i.  e. 


u.  (a.  +  u,  1 

i'  j  t- 


V.  u.  +  u.  u, 
II  I  k 
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(v)  There  exists  an  identity  element  Uq,  under  addition,  i.  e. 
u„  +  u.  =  u.  for  ciny  i  . 

0  j  j  ^  J 


(vi) 


There  exists  an  identity  element 


Uj,  Uj  =  Uj  for  any  j 


Uj^,  under  multiplication,!,  e. 


(vii)  For  each  element  u^  ttiere  exists  a  unique  inverse  with  respect 
to  addition,!,  e. 
u.  +  u.,  =  u_  . 

i  i‘  0 

(viii)  For  each  element  u.  )  there  exists  a  unique  inverse  with 

respect  to  multiplication,!,  c. 

u.  u. ,  =  u,  . 
i  i’  1 


The  finite  field  of  p  elements,  where  p  is  a  prime  number  may  be 
represented  by  u^  =  0,  Uj  =  1,  =  2  . . .  u^  ^  =  p-1,  in  which  addition 

and  multiplication  are  vhe  ordinary  arithmetic  operations  with  the  rule 
that  the  numbers  are  to  be  reduced  modulo  p . 

In  general,  a  Galois  field  of  p”^  elements  is  obtained  as  follows: 

Let  P{x)  be  a  given  polynomial  in  x  of  degree  m  with  integral 
coefficients;  and  let  F(x)  be  any  polynomial  in  x  v.'i*'h  integral 
coefficients.  Then  F(x)  may  be  expressed  as 

F(x)  =  f(x)  +  p.  q(x)  +  P{x)  Q(x) 

here  q(x)  and  Q{x)  may  be  any  polynomial  in  x  with  integral 
coefficients  and 

rt  \  .  .  2  ,  ,  m-1 

ffx)  =  a_  +  a,  x  +  a,,x  +  ...+a  ,x 
''01  2  m-1 

and  the  coefficients  a^,  a^,  ....  a^  j  belong  to  the  set  0,  1,2,  ...,p-l. 
This  relationship  may  be  written  as 
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F{x)  =  £(x)  mod  p,  P(x) 


and  f(x)  is  said  to  be  the  residue  of  F{x)  modulis  p  and  P(x).  It  p 
and  P(x)  are  kept  fixed  then  the  f(x)>s  form  p™  classes  of  functions. 
It  may  be  readily  verified  that  when  p  is  a  prime  number  and  P(x)  is 
irreducible  modulo  p.  that  is,  P(x)  cf nnot  be  expressed  in  the  form 

P(x)  =  Pj(x)  P^Cx)  +  p.  P3(x) 

then  the  classes  defined  by  the  f{x)’s  make  up  a  field. 

The  finite  field  formed  by  the  p”'  classes  of  residues  is  called  a 
Galois  field  of  order  p”^  and  is  denoted  by  GF(p”^).  The  p*”  classes 
are  the  same,  regardless  of  the  choice  of  Pfx),  subject  to  the 
restrictions  imposed  above,  and  the  field  GF(p”'),  always  exists  if  p 
is  a  prime  and  m  a  positive  integer.  The  classes  of  residues  can  be 
represented  by  the  different  possible  functions  f(x)  and  may  also  be 
denoted  by  Uq,  Uj,  U2,  ...»  u^_j  where  s  =■  . 

2 

To  illustrate,  we  shall  obtain  the  Galois  field  of  3  elements.  An 
irreducible  polyno  .nial  modulo  3  is  P(x)  =  1  +  x^.  Now  consider  the 
possible  functions  f(x).  These  are  of  the  form  a_  +  a,  x  where  a„  and 

U  i  w 

a^  are  elements  of  the  set  0,  1  and  Z.  Hence  the  elements  of  the  field 
are:  Uq  =  0,  Uj  =  1,  u^  =  Z,  =  x,  =  Zx,  =  1  +  x,  =  1  +  Zx, 

u.^  =  2  +  X,  Ug  =  2  +  2x.  There  is  a  further  theorem  that  all  the  elements 
or  marks  of  the  field  except  the  zero  element  Ug  can  be  represented  as 
the  powers  of  an  element  known  as  a  primitive  mark.  It  is  readily 
verified  that  y  =  1  +  x  is  a  primitive  mark.  For 
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y  =  1  +  X 

y  =  l  +  2xi-x‘'=2x  (since  P(x)  =  1  +  x'^  =  0) 

=  2x  +  2x^  =  2x  +  1  +  2(1  +  x^)  =  1  +  2x 

4  2  2  2 

y  =  4x‘^  =  x^  =  2  f  1  +  x'^  =  2 

y^  =  2  +  2x 

y^  =  2+  2x  2x+  2x^=.  4x+  2(1  +  x^)  =  4x=  X 

y^  =  X  +  x^  =  X  +  2  +  (1  +  x^)  =  2  +  X 

y  =  4  t;  1 

Both  the  representations  of  the  elements  of  the  field  are  important  in 
that  the  representation  in  terms  of  x  is  used  for  addition  and  the 
representation  in  terms  of  y  for  multiplication. 

An  irreducible  polynomial  P(x),  a  primitive  mark  and  the  addition 
and  multiplication  tables  for  GF(2^),  GF(2^)  and  GF(3^)  are  now 


presented. 


TABLE  6 


THE  GALOIS  FIELD,  GF(2^) 


2 

P(x)  =  1  +  X  •(-  X  ,  Primitive  mark  =  x . 
Addition  Multiplication 


0  12  3 


0  12  3 


0 


0  0  0  0 


1  1  Z  3 


2 

3 


3  1 


2 


TABLE  7 

THE  GALOIS  FIELD,  GF(2^  ) 

2  3 

P(x)  =  1  +  X  +  X  ,  Primitive  mark  =  x 
Addition  Multiplication 


0 

1 

2 

3 

4 

5 

6 
7 


0  1  2  3  4  5  ^7 

0  1  2  3  4  5  6  7 

0  3  2  5  4  7  6 

0  1  6  7  4  5 

0  7  6  5  4 

0  12  3 

0  3  2 

0  1 
0 


0 

1 

2 

3 

4 

5 

6 
7 


0  1  2  3  4  5  6  7 

QOOOOOOO 
1  2  3  4  5  6  7 

4  6  5  7  1  3 

5  12  7  4 

7  3  2  6 

6  4  1 

3  5 

2 
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TABLE  8 


THE  GALOIS  FIELD,  GF(3^) 


P(x) 

1  + 

2 

X  , 

Primitive  mark 

1  + 

X 

Addition 

Multiplicatio 

n 

0  1 

2 

3 

4 

5 

6 

7 

8 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

0  1 

2 

3 

4 

5 

6 

7 

8 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

2 

0 

4 

5 

3 

7 

8 

6 

1 

1 

2 

3 

4 

5 

6 

7 

8 

2 

1 

5 

3 

4 

8 

6 

7 

2 

1 

6 

8 

7 

3 

5 

4 

3 

6 

7 

8 

0 

1 

2 

3 

2 

5 

8 

1 

4 

7 

4 

8 

6 

1 

2 

0 

6 

1 

7 

2 

3 

r 

J 

7 

2 

0 

1 

3 

4 

6 

2 

6 

3 

4 

5 

6 

2 

8 

5 

7 
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3 

7 

3 

1 

8 

4 

8 

6 

The  use  of  these  fields  in  examining  the  s”  =  factorial  system 

is  ev:\ctly  analogous  to  the  use  of  0,  1,  2,  ,  p-1  for  the  p"  factorial 

system.  The  treatment  combinations  may  be  denoted  by  (Xj.x^, . .  .  ,x^) 
where  each  x^  can  take  one  of  the  values  0,  1,  2,  ....  s-1.  The 
(s”-l)/(s-l)  sets  of  (s-1)  degrees  of  freedom,  which  can  be  obtained  by 
partitioning  the  (s^-1)  degrees  of  freedom  into  main  effects  and 
interactions,  are  orthogonal,  and  the  responses  may  be  expressed  in 
terms  of  the  mean,  effects  and  interactions.  The  only  complication  is 
that  the  numbers  used  are  marks  of  the  Galois  field,  addition  and 
multiplication  being  defined  within  the  field. 
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C.  Confounding  in  Factorial  Experiments 


The  pcilormaace  of  a  comparative  experiment  requires  definition  of 
experimental  units  and  the  precision  of  conclusions  depends  on  the 
variation  among  the  units,  in  addition  to  other  things.  The  greater  the 
variation  among  units  the  higher  the  error  and  the  lower  the  precision. 

To  combat  this,  it  is  advantageous  to  group  the  units  into  what  are 
usually  called  blocks  of  units  and  to  design  the  experiment  so  that  only  the 
variation  among  units  within  blocks  enters  the  standard  error  of 
estimates.  The  smaller  the  block  size  the  more  uniform  the  units  in  the 
block  will  tend  to  be.  It  is  therefore  desirable  to  have  some  means  of 
reducing  the  size  of  the  block,  i.  e.-  the  number  of  units  in  each  block, 
and  thus  increase  precision.  For  this  purpose  the  device  of  confounding 
has  been  found  very  useful. 

Consider  a  simple  situation  of  three  factors.  A,  B  and  C  each  at 
two  levels,  the  eff  cts  and  interactions  being  defined  as  in  Table  V 
(apart  from  the  conventional  numerical  divisor). 

Suppose  that  the  eight  treatment  combinations  are  arranged  in  two 
blocks  according  to  their  sign  in  the  ABC  interaction.  The  two  blocks 
would  then  contain  the  following  trea<-ment  combinations: 


Block  1 

Block  2 

(1) 

a 

ab 

b 

ac 

Q 

be 

abc 
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TABLE  9 


EFFECTS  AND  INTERACTIONS  OF  THE  7?  EXPERIMENT 


(1) 

a 

b 

ab 

c 

ac 

be 

abc 

A 

-  i 

1 

-1 

1 

-1 

1 

-1 

1 

B 

-1 

-1 

1 

1 

-1 

-1 

1 

1 

AB 

1 

-1 

-1 

1 

1 

-1 

-1 

1 

C 

1 

-1 

-1 

-1 

1 

1 

1 

1 

AC 

1 

-1 

1 

-1 

-1 

1 

-1 

1 

BC 

1 

1 

-1 

-1 

-1 

-1 

1 

1 

ABC 

-1 

1 

1 

-1 

1 

-1 

-1 

1 

The  quantity  used  to  estimate  A  is  orthogonal  to  blocks  in  that  it  is 
given  by  ^(  -  (1)  +  a  -  b  +  ab  -  c  +  ac  -  be  +  abc)  and  of  the  four  treat¬ 
ment  combinations  entering  the  estimate  positively  two  ^re  in  each  block, 
and  likewise  for  the  four  treatment  combinations  entering  negatively. 

The  estimate  will  then  contain  none  of  the  additive  block  effects.  The 
same  is  true  of  the  other  main-effects  and  the  two-factor  interactions. 

The  three-factor  interaction  is  estimated  by 
^(  -  (1^  +  a  +  b  -  ab  +  c  -  ac  -  be  +  abc)  and  this  estimate  measures  not 
only  the  true  ABC  interaction  but  also  the  block  difference  (block  2 
minus  block  1).  It  is  not  possible  to  separate  the  true  interaction  from 
the  block  difference  and  the  interaction  and  block  difference  are  said  to 
be  completely  confounded  with  each  other.  Thus  the  three-factor 
interaction  cannot  be  estimated.  In  many  situations  it  is  known  that  the 
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high  order  interactions  are  trivial  and  therefore  can  be  used  as  blocking 
factors. 

The  set  of  treatment  combinations  in  the  block  of  a  confounded 
experiment  which  includes  the  control  treatment  is  called  the  intrablock 
subgroup. 

If  none  of  the  interactions  can  be  considered  trivial  and  smaller 
blocks  are  desired,  the  experiment  can  be  replicated  several  times  with 
a  different  effect  or  interaction  confounded  with  blocks  in  each  replicate. 

3 

For  example,  in  the  2  experiment  we  might  replicate  the  experiment 
four  times  confounding  ABC  with  blocks  in  the  first  replicate,  AB  with 
blocks  in  the  second  replicate,  AC  with  blocks  in  the  third  replicate 
and  BC  with  blocks  in  the  fourth  replicate.  Thus  each  interaction  may 
be  estimated  in  the  three  replicates  in  which  it  is  unconfounded.  This 
type  of  confounding  is  known  as  partial  confounding. 

The  rule  of  the  generalized  interaction  for  2”  experiments  is  that  if 
effects  or  interactions  represented  by  X  and  Y  are  confounded,  then  so 
is  XY,  obtained  by  multiplying  the  symbols  together  equating  any  letter 
which  is  squared  to  unity.  The  rule  of  the  generalized  interaction  for  the 
3”  system  is  that  if  pairs  of  degrees  of  freedom  corresponding  to  X  and 
Y  are  completely  confounded,  then  so  are  the  pairs  of  degrees  of  freedom 
corresponding  to  XY  and  XY^  where  any  letter  cubed  is  equated  to 
unity.  By  adopting  the  rule  that  in  any  symbol  the  power  of  the  first 
letter  should  be  unity  a  complete  and  unique  specification  of  all  effects 
and  interactions  is  achieved. 

An  example  of  the  use  of  this  symbolism  is  the  following.  Suppose  a 
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3  experiment  is  to  be  arranged  in  blocks  of  three.  In  any  system  of 

confounding  four  pairs  of  degrees  of  freedom  must  be  confounded  with 

2  2 

blocks.  For  example,  if  AB  C  and  AC  are  confounded,  so  is 
AB^  C  X  AC^  =  A^  B^  =  A^  B^  =  A^  b'^  =  AB 

7  7  7  %  7  ^  7  7  44 

and  AB'^C  x  AC  KC  =  ABC  =  B  C  =  B  =  BC  . 

The  composition  of  the  blocks  is  easily  obtained  from  the  definition 
of  the  effects  and  interactir  -.s.  In  th.?  above  example  there  are  nine 
blocks  given  by  the  solutions  of  the  equations 

Xj  +  2x2  +  x^  =  i  (mod  3^ 

Xj  +  2x2  =  j  (mod  3) 

where  i  and  j  each  take  on  the  values  0,  1  and  2. 

The  rule  of  the  generalized  interaction  for  the  experiment  is  that, 
if  effects  or  interactions  denoted  by  X  and  Y  are  completely  confounded 
with  blocks,  then  so  are  the  (p-1)  sets  of  (p-1)  degrees  of  freedom 
denoted  by  XY,  XY^,  ...»  XY^  where  any  letter  raised  to  the 
powe-  is  to  be  replaced  by  unity  and  the  resultant  symbol  is  to  be  raised 
to  such  power  as  makes  the  first  letter  in  it  have  a  power  of  unity.  This 
may  be  proved  as  follows:  Let  X  correspond  to  the  equations 

+  02*2  ■*■•••  “n’^n  “  ~  ’’  ”  ‘  ‘ (mod  p) 

and  y  to  the  equations 

Pi  Xj  +  ^2^2  +  •  •  •  +  PjjXjj  =  0,  =  1,  =  .  .  .  =  (p-1)  (mod  p)  . 

Because  X  and  Y  are  confounded  completely  with  blocks,  the  treatment 
combinations  of  any  one  block  satisfy  the  equations 
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ajXi+a2X2+...  ^  =  i  (mod  p) 

P1X1+P2X2+...  =  j  (modp) 

where  i  and  j  are  each  one  of  the  numbers  0,  1,  (p-1)-  For 

these  treatment  combinations  the  equations  may  be  combined  to  give 

(gj  +  XPi)Xj  +  (a^  +  ‘  ^^n^^n  =  ^  P) 

where  X  can  take  on  any  value  from  1  to  (p-1)  and  the  coefficients  on 
both  sides  of  the  equation  must  be  reduced  modulo  p.  This  equation 
corresponds  to  the  symbol  XY^ .  Thus,  the  treatment  combinations  of 
any  block  take  on  a  constant  value  for  any  one  if  the  equations  corre¬ 
sponding  to  XY^  where  X  is  any  value  from  1  to  (p-I).  The  effect  or 
interaction  XY^  is  therefore  confounded  with  blocks  for  these  values  of 
X. 


D.  Fractional  Replication 

A  complete  factorial  experiment  investigating  all  possible  combinations 
of  all  the  ^Levels  of  the  different  factors  will  invol’^e  a  large  number  of 
trials  when  the  number  of  factors  is  five  or  more.  When  the  number  of 
factors  is  large  the  number  of  trials  required  may  even  become 
prohibitive.  One  is  therefore  led  to  consider  the  economy  of  space  and 
material  which  will  be  attained  by  using  only  a  fraction  of  the  possible 
number  of  treatment  combinations  at  the  expense  of  losing  some 
information  inherent  in  a  complete  replicate.  The  general  process  by 
which  information  can  be  obtained  from  less  than  a  full  replicate  of  a 
factorial  experiment  is  known  as  fractional  replication. 

Suppose  that  three  factors,  A,  B  and  C,  each  having  two  levels  are 
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under  investigation  and  it  is  known  that  these  factors  do  not  interact. 

The  relation  between  the  true  responses  and  the  effects  and  interactions 
can  be  presented  in  tabular  form  as  follows; 


TABLE  10 

RELATION  BETWEEN  RESPONSES  AND  EFFECTS  AND  INTERACTIONS 

IN  A  EXPERIMENT 


(A 

1 

Iab 

Iac 

Ibc 

Iabc 

(1) 

+ 

- 

- 

+ 

- 

+ 

+ 

- 

Zl 

+ 

+ 

- 

- 

- 

- 

+ 

+ 

b 

+ 

- 

+ 

- 

- 

+ 

- 

+ 

ab 

+ 

+ 

+ 

+ 

- 

- 

- 

- 

c 

- 

- 

+ 

+ 

- 

- 

+ 

ac 

+ 

+ 

- 

- 

+ 

+ 

- 

- 

be 

+ 

- 

+ 

- 

+ 

- 

+ 

- 

abc 

+ 

+ 

4* 

+ 

+ 

T 

+ 

+ 

Suppose  that  only  the  four  treatmeiit  combinations  that  enter  the  ABC 
interaction  negatively  are  considered;  namely  (1),  ab,  ac  and  be.  It 
is  clear  from  the  table  that  it  is  impossible  to  separate  the  mean  p,  from 
the  ABC  interaction.  Similarly  the  A  effect  cannot  be  separated  from 
the  BC  interaction,  the  B  effect  ca  mot  be  separated  from  the  AC 
interaction  and  the  C  effect  cannot  be  separated  from  the  AB  inter¬ 
action.  The  estimating  equations  for  this  plan  are; 
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lA  -  jAEC 

I  (A  -  BC) 
i  (B  -  AC) 
i  (C  -  AB) 


^  £*(1)  +  ab  +  ac  +  be  J 

^  /f  -  (1)  +  ab  +  ac.  -  be  J 

^  T-  (1)  +  ab  -  ac  +  be  J 

^  r  -  ^1)  -  ab  +  ae  +  be  7 


With  only  the  four  trials,  A  is  completely  eonfounded  or  aliased 
with  EC,  B  with  AC,  with  AB  and  p  with  ABC. 

ff  the  factors  do  not  interact  all  the  interactions  may  be  neglected 
and  the  estimating  equations  can  then  be  used  to  estimate  the  mean  and 
the  three  main  effects,  and  these  estimates  are  uncorrelated. 

The  treatments  chosen  for  the  1/2  replicate  were  selected  as  those 
which  entered  the  ABC  interaction  negatively.  The  selection  could  have 
been  those  which  entered  the  ABC  interaction  positively.  For  most 
purposes,  it  does  not  matter  which  of  the  two  halves  of  the  experiment  is 
chosen. 

If  by  convention,  p  is  denoted  by  I,  the  confounding  relation  may  be 
expressed  as 

I  =  ABC  . 

This  relation  is  known  as  the  defining  contrast  or  identity  relationship 
where  the  equal  sign  is  used  to  denote  "completely  confounded  with". 

The  remaining  three  confounding  relations  may  be  written  as 


A  =  BC 
B  =  AC 
C  =  AB 
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These  relations  may  be  obtained  from  the  identity  relationship  by 
multiplying  both  sides  by  an  effect  of  interest  with  the  rule  that  any 
letter  which  is  squared  is  to  be  replaced  by  unity.  Thus  I  =  ABC  when 
multiplied  by  A  gives 

A  =  .f  ^BC  =  BC  . 

3 

It  should  be  noted  that  the  1/2  replicate  of  the  2  experiment  given 

by  the  identity  relationship  I  =  ABC  consists  of  the  same  treatment 

3 

combinations  as  one  of  the  blocks  of  a  2  experiment  in  two  blocks  of 
four  plots  each  where  ABC  is  confounded  with  blocks. 

The  identity  relationship  of  a  quarter  replicate  of  a  2”  experiment 
is  of  the  form 

I  =  X  =  Y  =  XY 

where  X,  Y  and  XY  are  higher  order  interactions  and  XY  is  the 
generalized  im'-raction  of  X  and  Y. 

If,  for  a  1/3  replicate  of  the  3^  experiment,  the  identity  relation¬ 
ship  is  given  by 

I  =  ABC 

the  following  confounding  relations  may  be  generated: 

A  =  AB^C^  =  BC 

B  =  AB^C  =  AC 

C  =  ABC^  =  AB 

2  2  7 

AB  =  AC  =  BC*^ 
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In  general  a  1/p”^  replicate  of  a  factorial  experiment  may  be 
specified  by  an  identity  relationship  of  the  form 

I  ^  X  Y  =  XY  =  XY^  =..,=.  XY^"'  =  Z^XZ  =  XZ^=...  =  XZ^"^ 

=  YZ  =  YZ^  =  . . .  =  YZ^'^  =  XY^Z®  Tr,  s  =  1,  2,  .  ,  (p-l)7 

=  etc. 

v-fhf^vp  there  are  t  independent  contrasts  X,  Y,  Z,  etc. 

For  a  more  detailed  discussion  of  factorial  experiments  one  may 
refer  to  the  texts  by  Kempthorne  (1952)  and  Davies  {1954}. 
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IV.  ORTHOGONAL  MAIN-EFFECT  PLANS 


The  experimental  plans  which  are  developed  in  this  chapter  and 
subsequently  presented  in  the  catalogue  are  called  orthogonal  main-effect 
plans  as  they  permit  the  estimation  of  all  main  effects  without  correlation, 
when  all  interactions  are  negligible. 

The  most  commonly  used  factorial  experiments  involve  factors  which 

all  occur  at  the  same  number  of  levels.  These  experiments  are  known 

as  symmetrical  factorial  experiments.  A  good  deal  is  already  known 

about  the  construction  of  orthogonal  main-effect  plans  for  symmetrical 

factorial  experiments  although  a  comprehensive  catalogue  of  such  plans 

has  never  been  published.  There  are  a  great  many  experimental 

situations  which  involve  factors  that  do  not  all  occur  at  the  same  number 

of  levels.  These  experiments  are  known  as  asymmetrical  factorial 

experiments.  Heretofore  the  standard  technique  for  constructing 

orthogonal  main-effect  plans  for  asymmetrical  factorial  experiments  has 

been  to  combine  two  or  more  orthogonal  main-effect  plans  for  different 

symmetrical  experiments.  Hence  in  order  to  construct  an  orthogonal 

4  3 

main-effect  plan  for  the  3  x  2  experiment  one  would  combine  the  plan 

4  3 

for  the  3  experiment  in  nine  trials  with  the  plan  for  the  2  experi¬ 
ment  in  four  trials  to  obtain  the  required  plan  in  thirty- six  trials.  This 
^rocedure  often  requires  more  trials  than  the  experimenter  can  afford  to 
make. 

The  orthogonal  main-effect  plans  developed  in  this  chapter  for  both 
symmetrical  and  a  .symmetrical  factorial  experiments  require  the  least 
number  of  trials  that  has  yet  been  attained  for  such  plans.  For  meny 
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experiments  the  suggested  plans  are  so  highly  fractionated  that  there  •'re 
few  if  any  degrees  of  freedom  available  for  the  estimation  of  experimental 
error.  In  such  situations  one  must  use  an  estimate  of  experimental  error 
which  is  (i)  known  from  previous  experience,  (ii)  derived  from  some  of 
the  degrees  of  freedom  available  for  estimating  main-effects  which  prior 
knowledge  indicates  are  negligible  or  (iii)  appro:timatv,d  by  a  procedure 
which  utilises  a  graph  and  is  known  as  the  half-normal  plot  technique  of 
interpreting  factorial  experiments.. 

The  plans  consist  of  th  reatment  combinations  which  permit 
uncorrelated  main  effect  estimates.  The  treatment  combinations  are 
denoted  by  the  level  at  which  each  factor  occurs.  Thus  the  treatment 
combination  0112  in  an  experiment  on  four  factors  is  that  combination 
for  which  the  first  factor  occurs  at  its  first  level,  the  second  and  third 
factors  occur  at  their  second  levels  and  the  third  factor  occurs  at  cs 
third  level.. 


A.  Weighing  Plans 

The  problem  of  estimating  the  weights  of  small  objects  placed  on  a 
balance  scale  was  first  considered  by  Yates  (1935).  The  weighing 
problem  is  concerned  with  the  development  of  plans  for  estimating  the 
effects  of  two-level  factors  with  as  few  trials  as  possible.  Since  it  can 
je  assumed  that  the  weight  of  a  set  of  objects  is  the  sum  of  the  weights  of 
the  individual  objects,  all  interactions  may  be  presumed  to  be  absent. 
Hotelling  (1944)  constructed  optimum  (in  the  sense  of  minimum  variance) 
plans  for  estimating  the  weights  of  (N-1)  objects  with  N  weighings  on 
a  chemical  balance  scale.  He  proved  that  a  necessary  and  sufficif-nt 
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condition  for  attaining  an  optimum  weighing  plan  is  that  the  design 
matrix,  X  say,  be  a  Hadamard  matrix,  which  is  a  matrix  consisting  of 
I’s  and  -i’s  such  that  X*X  =  diagonal  (N,  N, ...,  N)  where  N  is  the 
number  of  weighings.  Paley  (1933)  proved  that  a  siificient  condition  that 
a  Hadama  1  Cl  liix  of  si-^e  N  exist  is  N  =  0  {mod  4),  with  the  exception 
of  N  =  2  which  is  a  trivial  case. 

Placket*^  and  Burman  (1946)  provided  what  is  effectively  a  complete 
solution  of  the  weighing  v  oblem  when  the  estimates  of  the  weights  are 
required  to  be  uucorrelated.  Most  of  the  plans  which  they  developed  can 
be  generated  by  a  cyclic  shifting  of  the  elements  of  one  treatment 
combination  successively  (N-ZjT  times  and  then  adding  the  control  treat¬ 
ment.  When  the  number  of  trials  N  =  0  (mod  4)  is  not  of  the  form 
N  =  2”  the  orthogonal  main-effect  plans  given  in  the  catalogue  have  been 
generated  by  cyclically  shifting  the  elements  of  the  treatment  combi¬ 
nations  presented  by  Plackett  and  Burman. 

j.  .  Plans  for  Symmetrical  Factorial  Experiments 

Orthogonal  main-effect  plans  can  be  constructed  easily  for 
symmetrical  factorial  experiments  involving  (s*^-l)/(s-l)  factors,  each 
having  s  levels,  with  s*'  treatment  combinations,  where  s  =  p*”  and 
p  is  a  prime  number.  The  (s^-l)/(s-l)  factors  can  be  represented  by 
n  factors  each  having  s(=  p*^)  levels  and  all  their  generalized  inter¬ 
actions.  Hence  one  need  only  choose  the  treatment  combinations  from  a 
complete  s”  factorial  plan  and  assign  one  of  the  (s*'- l)/(s- 1)  factors 
to  each  of  the  factors  and  interactions  of  the  s"  plan. 

Let  the  n  factors  of  the  s*'  factorial  plan  be  represented  by 
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Xp  X^,  •  •  •  ,  X^  and  their  generalized  interactions  by 

k.  X,  +  k_  Xt +.  .  .  +  k  X  where  k.  (i  =  i ,  2, . .  .  ,  n)  can  take  on  any  value 

of  the  Galois  field  GF(p”')  and  it  is  understood  that  the  coefficient  of 

the  first  factor  appearing  in  an  interaction  is  unity.  The  notation 

adopted  here  for  the  generalized  interactions  differs  from  the  standard 

notation  for  interactions  as  given,  for  example  by  Kempthorne  (1952),  in 

order  to  facilitate  the  presentation  which  follows  later., 

The  procedure  for  constructing  orthogonal  main-effect  plans  will  be 

illustrated  with  a  plan  for  four  factors  A,  B,  C  and  D,  each  having 

three  levels  with  nine  treatment  combinations.  In  this  example  s  =  3, 

n  =  2  and  {s*^-l)/{s-l)  =  4.  The  four  factors  can  be  represented  by  two 

2 

factors  Xj  and  of  the  3  factorial  experiment  and  their  general¬ 
ized  interactions  Xj  +  X^  and  Xj  +  ^X^-  The  treatment  combinations 
which  comprise  the  orthogonal  main-effect  plan  are 

0  0  0  0 
0  112 
0  2  2  1 
10  11 
112  0 
12  0  2 
2  0  2  2 
2  10  1 
2  2  10 

The  interactions  which  are  members  of  the  defining  contrast  (identity 
relationship)  may  be  determined  by  choosing  those  interactions  whose  X 
representation  equals  0  (mod  3).  The  generators  of  the  interactions  in 
defining  contrast  for  the  example  given  above  are  ABC^  and  ACD, 
since  the  X  representation  of  ABC^  is  +  X2  +  2(Xj  +  X^)  =  0 
(mod  3)  and  of  ACD  =  Xj  +  (Xj  +  X2)  +  (Xj  +  ZX^)  =  0  (mod  3). 
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Some  plans  which  may  be  constructed  by  this  method  are  given  in 


Table  H. 


TABLiE  il 

INDEX  OF  SOME  MAIN-EFFECT  PLANS 


Number 
of  levels 

Maximum  Number 
of  factors 

Number  of 
observations 

2 

3 

4 

2 

7 

8 

2 

15 

16 

2 

31 

32 

2 

63 

64 

3 

4 

9 

3 

13 

27 

3 

40 

81 

A 

5 

16 

4 

21 

64 

5 

6 

25 

7 

8 

49 

8 

9 

64 

9 

10 

81 

The  orthogonal  main-elfect  plans  with  treatment  combinations 
which  accommodate  up  to  (s*^-l)/(8-l)  factors  can  be  augmented  to 
yield  orthogonal  main-effect  plans  with  2  s*'  treatment  combinations. 
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The  augmented  plans  can  accommodate  up  to  C  2  (s^- l)/(s- 1)  -  I J 
factors,  each  having  s  =;  p^  levels.-  In  order  to  illustrate  the  theory- 
underlying  the  augnie’-.tation  procedure  some  preliminary  lemmas  are 
now  developed. 

Let  u„,u,, .  . .  ,u  ,  represent  the  elements  of  the  Galois  field 

U  1  3—1 

ni  Z  Z  Z 

GF(p  )  and  let  Uq,  u^,  ...,u^  ^  represent  the  squares  of  the  elements 

of  GF{p”^).  The  set  of  squared  elements  of  GF(p”^)  will  be  denoted  by 

GF^(p*^).  It  is  easily  ve-  Ified  that  apart  from  the  0  element,  the  set 
Z  m 

GF  (p  )  forms  a  cyclic  Abelian  group  under  multiplication.  It  follows 

Z  tti 

from  the  cyclic  property  that  (i)  when  p  =  2,  GF  (p  )  contains  each  of 
the  elements  of  GF{p”^)  and  (ii)  when  p  is  an  odd  prime,  the  elements 
of  GF^tP*^)  comprise  a  subset  of  (s  +  l)/2  distinct  elements  of 
GF(p^),  where  one  element  occurs  once  and  (s-l)/2  elements  are 
duplicated.- 

Consider  one  of  the  factors  in  a  main-effect  plan  in  which  each 

n—  1  n 

X.  has  s  levels  each  occurring  s  ”  times  in  a  total  of  s  treatment 

2 

combinations.  Let  X^'  be  a  pseudo-factor  obtained  by  squaring  the 
levels  of  X..  The  following  lemmas  can  now  be  presented: 

2 

Lemma  1:,  When  p  is  an  odd  prime,  X^  +  kXj^  (k  an  element  of 
GF(p™) )  contains  (s  +  l)/2  distinct  levels,  one  level  occurring  s*^ 
times  and  {s-l)/2  levels  occurring  2  s”  ^  times  in  s*^  treatment 
combinations. 

Z  n  1 

Lemma  2:  When  p  =  2,  X.  contains  e’  ..  of  the  s  levels  s  times. 

Lemma  3:  When  p  =  2,  X?  +  kX^  IJc  ^  Oj  contains  s/2  distinct  levels 
each  occurring  2  s^  ^  times. 


64 


Lemma  3  can  be  proved  as  tcllows.  Let  range  over  the  elements  of 
GF(p*^)  which  represent  the  s  levels  of  X^^.,  As  ranges  over  the 
elements  of  the  field  so  does  x^^  +  k  where  k  is  an  element  of  GF(p^). 
Also  if  +  k  =  Xj  (mod  ?.)  then  x^  +  k  =  Xj^  (mod  2).  Hence 
x.(x.  +  k)  =  X.  X.  and  x.(x.  +  k)  =  x.  x..  Thus  whatever  values  of 

i  11  j  '  i  j 

Xi(Xi  +  k)  are  achieved  they  are  achieved  for  at  least  two  values  of  x^. 

It  will  now  be  shewn  that  the  values  of  x.(x^  +  k)  are  achieved  for 
exactly  two  values  of  x..  I  et  y  be  the  generator  of  the  field  and  let 
x^  =:  y°’  and  k  =  y^.  Thus  +  k)  =  y**(y^  +  y^). 

Suppose  that 

y“(y“  +  y^)  =  y''(y’'  +  y'*) 

where 


a  a  y 
y  4  r 


and 


a  p 

y  +  y 


4y^ 


Hence 

(y*^)^  +  y“y^  =  (y^)^  +  y^y^ 

(y*^  +  y^l^  +  (y*^  +  y^ )  y^  =  o 

(y®  +  y'" )  (y®  +  y^  +  y^ )  =  0  . 

This  implies  that  either  y®  +  y^  =  0  and  therefore  y**'  =  y^  which  is  a 
contradiction  or  that  y®  +  y^  +  y^  =  0  and  therefore  y**^  +  y^  =  'p'  which 
»s  a  contradiction.  Hence  the  values  of  x.(x.  +  k)  are  achieved  for 
exactly  two  values  of  x.  and  Lemma  3  is  proved. 

2 

Lemma  4:'  The  factor  represented  by  X.  +  k.  X.  +  S  k.X.  ,  where  at 

■  t  t  1  -j-  J  J 

least  one  k.  4  0,  contains  each  of  the  s  levels  s*^  ^  times. 
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Lemma  5;  The  levels  of  X.  +  k,  X.  +  kTX.  which  occur  in  a  plan  with 

1  1  i  2  j 

the  level  of  a^X^  where  k^,  k^,  a^  and  are  elements 

of  GF(p*^)  and  a^  ^  0  are  given  by  the  values  of 
2 

X.  +  k,  X.  +  k^x.  +  cfa,  x. -fa^x. }  -  cu^  where  k-,  +  caT  =  0  and  x. 
ili2j'li2j'  t  2  2  1 

ranges  over  the  elements  of  GF(p  ). 

Proof:  When  a,  X.  +  a^X.  takes  on  the  level  then  a,  x.  +  a_x.  =  u. 

Ii2j  t  li2jt 

and  thus 


Xj  =  (u^  -  ajx.  )/a 


v-2. 


•2' 


Hence  the  levels  of  the  factor  X.  +  k,  X.  +  k_X.  which  occur  with  the 

1  1  1  2  j 

lev  -I  u^  of  aj^Xj^  +  a^Xj  can  be  represented  by 

xf  +  k^  X.  +  k^x.  =  xf  +  kj  X.  +  k^Cu^  -  a  J  x.)/a2 


=  x^  +  (kj  -  k2aj/a2)x.  +  (k^/a2)u^  . 


Since  k^  +  ca^  =  0.  then  c  =  -k^/a^  • 

Thus 

2  2 

X.  +  (k,  -  k,  a,  /a_)  X.  +  fk_/a,)  u  =  x.  +  k,  x.  +  k_x.  +  c(a,  x.  +  a,x.)  -  cu.  , 

1  ‘  1  2  1  2'  i  •  2  2'  t  1  1  1  2  J  •  1  1  2  y  t’ 

and  the  lemma  is  proved. 

Two  factors  X.  and  X^  are  said  to  be  orthogonal  to  each  other  if 
each  level  of  X.  occurs  the  same  number  of  times  with  every  level  of 
X- .  Two  factors  Xj^  and  Xj  are  said  to  be  semi-orthogonal  to  each 
other  if  fi)  for  p  an  odd  prime,  one  level  of  X^  occurs  s*^"^  times 
and  (s-l)/2  levels  of  X.  each  occur  2s”'^  times  with  each  level  of 

J 

■J 

X.  and  (ii)  for  p  -  2,  b/2  levels  of  X^  each  occur  2  3*^"  times  with 

each  level  of  X.  . 

1 
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It  follows  from  Liemmas  1,  3,  and  5  that  when  p  is  an  odd  prime 


or  when  k,  -  k^a./a^  4  G,  tlicxi  a.  X.  r  a,X.  is  semi-orthogonal  to 
2 

X.  4-  k,  X.  +  k^  X.  .  It  follows  from  Lemmas  2  and  5  that  when  p  =  2 

1  1  1  2  j  ^ 

and  kj^  -  k^aj/a^  =  0  then  +  a^X.  is  orthogonal  to 

2 

+  kj  Xj^  +  k2Xj  .  Employing  an  argument  similar  to  that  used  in 

2  2 

Lemma  5  it  can  be  deduced  that  kX.  +  k,  X.  +  X.  and  kX.  +  X.  +  X. 

i  1  1  j  1  2  1  1 

are  orthogonal  to  each  other  when  k^  ^  k^- 

Lemma  5  can  be  generalized  to  include  more  than  two  factors  as 
stated  in  Lemma  5a. 

2 

Lemma  5a:  The  levels  of  X.  1  k.  X.  +  S  k.X.  which  occur  in  a  plan  with 

^  I  j^i  J  J 

the  w  level  of  a.  X.  +  S  a.  X.  are  given  by  the  values  of 

^  "  j^i  J  J 


X.  1  k.  X.  +  2  k.  X.  +  c(a.  x.  +  S  a.  x. )  -  cu 

t  1  I  J  1  W  4  J  j'  t 


where  k.  +  ca.  =  0  for  all  j  ^  i.  If  the  a.  and  the  k.  are  not  ol  such  a 
J  J  J  J 

form  that  k^  +  ca^  =  0  for  all  j  ^  i  and  some  c  contained  in  GF(p”') 
then  the  two  factors  are  orthogonal. 

Lemma  6:  When  p  is  a  prime  the  complements  in  GF(p”*)  to  the 
elements  in  GF^^p*^)  are  the  set  of  elements  in  GF^(p”')  each  multi¬ 
plied  by  an  element  of  GF(p”^)  which  is  not  an  element  of  GF“(p^). 

Z  m 

If  the  set  of  elements  in  GF  {p  |  id  their  set  of  complements  are 
taken  together  in  ore  set  the  elements  of  GF(p”^)  are  obtaine'^. 


Proof:  From  abstract  group  theory  (see  Birkhoff  and  MacLane  (1953) ) 
we  employ  a  lemma  which  states  that  two  right  cosets  of  a  subgroup  are 
either  identical  or  v-ithout  common  elements.  Now  the  elements  of 
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Z 

GF  (p“  )  form  an  Abolian  subgroup  of  the  elements  d'  GF(p  ).  Hence 

multiplying  each  element  oi  GF^(p^)  by  an  element  of  GF(p*^)  which 

is  not  an  element  of  GF  (p  )  yields  the  complementary  set  to  GF^{p^). 

It  is  clear  from  Lenima  2  that  when  p  =  2  the  set  complementary  to 
Z  m 

GF  (p  )  is  the  null  set.. 

We  can  now  present 

Theo-em  1:  There  exists  a  main-effect  plan  for  C  2(s^‘-l)/(s-l)  -ij 
factors,  each  at  s  =  p*^  levels,  with  2s''  treatment  combinations. 

Proof:  ’'n  order  to  facilitate  the  presentation  of  the  proof  of  Theorem  1, 

2 

let  n  =  2.  First  construct  an  orthogonal  main  effect  plan  for  (s  ~I)/(e-1) 
factors  each  at  s  levels  in  s^  trials,  represented  by  the  two  factors 
Xj  and  and  their  generalized  interactions  Xj  +  X^i  Xj  +  2X2* 

+  (s-l)X2-  To  these  add  [ [s  -l)/(s-l)-l7  factors  represented  by 
Xi  +  Xp,  f  Xj  X2,  Xi  +  2X1  +  X,,  . .  ,  Xj  1  (s-l)Xj  +  X2.  Those 
/* 2 (s*'- 1) /(s- U  -  1  7  factors  in  s”  observations  repreient  the  first  half 
of  the  main-effect  plan. 

Note  from  the  proceeding  lemmas  that  when  p  is  a  prime  number, 

X,  +  3'j^X2  and  X^  +  k^^Xj  +  X2  aVe  semi-orthogonaJ  and  also  that  X2 

2  TTi 

and  Xj  +  Kj^Xj  V  X2  are  semi-orthogonal  for  all  a.  and  k.  in  GF(p  ; 
except  a^^  =  0,  All  other  pairs  of  factors  are  clearly  orthogonal.  If 

p  =  2  and  (k.  -  a^/a.)  =  0,  then  a,  Xj  +  a.  X2  and  +  k.  X^  +  X2  are 
orthogor.al. 

The  second  half  of  the  plan  is  chosen  so  that  the  pairs  of  factors  which 
are  orthogonal  in  the  first  half  are  also  orthogonal  in  the  second  half  and 
pairs  oi  factors  wliich  are  semi-orthogonal  in  the  first  half  are  sei.ii- 
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orthogonal  in  a  complementary  manner  in  the  second  half.-  The  factors  in 
the  second  half  which  correspond  to  the  factors  of  the  first  half  can  be 


denoted  by 

XpX2,,Xj  +  X2+bj,  Xj  +  2X2+b2,  •  .  •  ,  Xj  +  (s- IjX^  +  bg_  j ,  kX^  +  X^, 

kXj  •}  Xj  +  X,  +  Cj,  kXj  H-k^Xj  f  X2+  c,, .  .  .  ,  k X ^  t  k^g_  ^  X,  +X2+  Cg_j 

where  the  coefficients  b  j,  b^  j  •  •  »  ^g_  1 »  >  ^2’  ■■'’^s-l’^l’'^2’‘''’*^s  1 

arc  to  be  determined. 


From  Lemma  5,  it  is  seen  that  the  levels  of  Xj  +  which  occur 


with 


the  level  of  X~  are  given  by  the  values  of  +  u^  where 


takes  on  the  values  of  the  elements  of  GF(p  ).  Without  loss  of 
generality  we  may  let  P  prime,  the  values 

of  kX^  +  X^j  where  k  is  an  element  of  GF{p*^)  but  not  an  element  of 
2  rti 

GF  (p  '),  which  occur  with  the  =  0  level  of  are  given  by  the 

2  2  2 
values  of  kx^  .  As  shown  in  Lemma  6,  kx^  complements  Xj,  . 

Thus,  when  p  is  an  odd  prime  k  can  take  on  the  value  of  any 

element  in  GF(p“  )  which  is  not  an  element  of  GF  (p  ).  If  p  =  2  it 

is  cleM  r  from  Lemma  2  that  k  =  i. 

A  method  for  determining  the  constants  bj,  b^, .  .  .  .b^  j. 


m 


kj,  k^,  ...jk^  p  Cj,,C2»---,c^  j,  when  s  =  p  and  p  is  an  odd 

prime  is  now  presented.  In  order  that  the  levels  of  kXj  +  X2  which 

Occur  with  the  0  level  of  X^  +  +  b^^  be  the  complements  of  the 

2 

levels  of  X,  +  X,  which  occur  with  the  0  levels  of  X,  +  a.  Xt,  b. 

12  1  i  2  i 

must  be  such  that  the  values  which  kx,  -  (1/a.  )x,  -  b./a.  takes  when 

1  '  1'  1  11 

Xj  ranges  over  the  field  GF{p*^)  complements  the  values  which 
2  2 

Xj  -  (l/a^)Xj  takes.  Now  Xj  -  (l/aj^)Xj  consists  of  one  element  of 
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GF(p”^)  occurring  once  and  {s-l)/2  elements  occurring  twice.  Let  the 
unique  element  of  GF(p  *)  be  u^.  Then  x“  -  (1/a. )  =  Uj  must  have 

only  one  solution  as  ranges  over  the  elements  of  GF(p*^).  Thus 

2  2  2 
1/a.  +  4u,  =  0  and  hence  4u,  =  -  1/a.  .  Since  kx,  -  (l/a.  )x,  -  b./a. 

il  1  1  i'l'lii 

must  complement  Xj  -  {l/a^)xj  the  equation 

kx2-(l/a.)xj-b./a.  =  uj 

must  also  have  only  one  solution. 

Therefore 

1/a?  +  4k{b./a.  +  u,  )  =  0. 

i  'll  1 

2 

Substituting  4Uj  =  -  1/a.  in  this  equation  and  solving  for  b.  we  get 

b.  =  (k-  l)/4ka..  (1) 

2 

To  find  the  levels  of  +  d^Xj  +  which  occur  with  the  0  levels 

of  X.,  note  that  there  exists  an  element  of  GF(p^),  u,  say,  such  that 
2 

x^  +  dj^Xj  =  u^  has  only  one  solution. 

2  2 

Thus  d.  +  4u_  =  0  and  hence  4u,  =  -  d.  .  In  order  that  the  levels  of 

1  2  2  1 

kX^  +  +  X^+c.  which  occur  with  the  0  levels  of  X^  complement 

2  2 

those  given  by  x^  +  d.Xj,  then  kxj  +  +  Cj^  =  u^  must  have  only  one 

2 

solution..  Substituting  4u,  =  -  d.  in  this  equation  and  solving  for  c.  we 

w  1  1 

get 

c^  =  k?/4k  -  d.^/4.  (2) 

2 

To  find  the  levels  of  Xj  +  d.  Xj  +  X2  which  occur  with  the  0  levels 

of  Xj  +  Sj^X^  note  that  there  exists  an  element  of  GF(p^),  u^  say, such 
2 

that  Xj  f  (d.  -  1/a.  )xj  =  u^  has  only  one  solution. 
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Thus 


(d.  -  l/a.)^+4u-  =  0  and  4u,  =  -  {d.  -  1/a.)^  . 

'  1  i  3  3  '  1  i' 

2 

Since  kxj^  +  (k^  -  J/a^)Xj  +  (c^  -  h.Ja.^  must  complement 
2 

Xj  +  (dj^  -  l/aj^)Xj  ,  the  equation 

kXj  ^^kj  -  l/a.)Xj  +  (Cj^  -  b./a^)  =  u^ 

must  also  have  only  one  solution  as  Xj  ranges  over  the  dements  of 
GFCp"^).  Therefore 

(k.  -  1/a.}^  ,  4k  [  (c.  -  b./a.)  -  U3]  =  0  . 

2 

Substituting  4  0^=  "  -  1/a^)  and  equations  (1)  and  (2)  into  this 

equation  we  get 

k.  =  k  d.  .  (3) 

Hence  equation  (2)  can  be  rewritten  as 

c.  =  d.^  (k-l)/4  .  (4) 

Thus  k  is  determined  by  choosing  an  element  of  GF(p*^)  which  is 
Z  m 

not  an  element  of  GF  (p  ).  By  letting  a.  =  1,  2,  . . . ,  s-1  we  can 
determine  bj,  b^,  • .  •  .  b^  ^  from  equation  ^).  Then  setting 
d;  =  1,  2,  s-l  we  determine  k^,  k^t  •  •  • »  j  from  equation  (3) 

and  Cj,  c,,  . . ,  Cg  j  from  equation  (4). 

The  procedure  employed  above  cinnot  be  applied  when  p  =  2  since 

Z  m 

Xj  +  cXj  consists  of  s/2  elements  of  GF(2  ),  each  occurring  twice. 

Thus  there  exists  no  element  u  such  that  Xj  +  cx^  =  u  must  have  only 

one  solution. 
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We  deduce  from  Lemma.  2  that  when  p  =  2,  then  k  =  1.  In  order 

that  the  levels  of  ’  X2  which  occur  with  the  0  levels  of 

Xj  *  a.  X,  +  b.  (a.  =  1,  2,  3;  . .  .  ,  s-1)  compieinent  the  levels  of 
1  1^1*1 

2 

X,  X-.  which  occur  with  the  0  levels  of  X,  +  a.  X~  then  tlic  levels 

1  ^  1  i  Z 

CT\ven  by  -  (l/a/;Xj  -  h./a^^  must  complement  the  Icvciti  given  by 

Z  m 

X.  -  {l/aj^)Xj  when  x.^  ranges  over  GF(2  ).  It  is  easily  verified  that 

can  be  any  one  of  the  ^  elements  of  GF{2’^)  which  are  not 

given  by  x“  -  (l/a^)Xj  .. 

2 

In  order  that  the  levels  of  X.  t  k.  X^  -r  X-,  +  c.  which  occur  with  the 

111  2  i 

0  levels  of  X^  complement  the  levels  of  Xj  +  d.  X^  +  which  occur  with 

the  0  levels  of  X^.  then  the  values  given  by  Xj  +  k.  x^  +  must  comple- 

2 

ment  the  values  given  by  x,  -i-  d.  x, .  It  can  be  shown  that  k.  =  d.  and  c. 

•'ll!  Ill 

can  be  any  one  of  the  2*^  ^  elements  of  GF(2*^)  which  are  not  given  by 

2 

the  values  of  x,  d.  x,  . 

1  i  1 

2 

By  finding  the  values  of  which  occur  with  the  0 

levels  of  Xj  +  ^  which  complement  the  valnos  of 

2 

X,  T  d.  X,  +  X,  that  occur  with  the  0  levels  of  X,  +  a.  X-.,  a  set  of  b. 

I  i  1  r  I  1  2  i 

and  which  satisfy  all  the  requirements  to  have  the  second  half  of  the  plan 

complement  the  first  half  of  the  plan  can  be  determined. 

When  n  >  2  the  same  procedures  will  yield  the  desired  plans  if 

L^mma  5u  is  utilized  in  place  of  Lemma  5.  Thus  the  theorem  is  proved. 

The  orthogonal  main-effect  plans  for  £"2(5” -l)/(s-l)  -  1  J  factors 

each  at  s  =  p”'  levels  with  2  s"  treatment  combinations  which  are 

7  25  9 

included  in  the  catalogue  are  the  following;'  3  in  18,  3  in  54,  4  in  32, 

II  7 

5  in  50.  Bose  and  Buah  (1952)  have  constructed  the  plans  for  3  in  18 


9  n 

and  4'  in  32  by  other  procedures  and  have  shown  that  ^72  {s  - l)/(s- 1)  -  1  J 

is  the  maximum  number  oi  factors,  each  at  s  levels,  that  can  be 

accommodated  in  an  orthc^  anal  main-effect  plan  with  2  s*^  treatment 

combinations. 


C.  Condition  of  Proportional  Frequencies 

In  the  complete  factorial  experiment  the  levels  of  a  factor  occur 
equally  frequently  with  each  of  the  levels  of  any  other  factor.  This  condi¬ 
tion  is  sufficient  to  allow  uncorrelated  estimates  of  all  effects  ".id  inter¬ 
actions.  This  condition  is  also  sufficient  to  allow  uncorrelated  estimates 
of  the  main  effects  in  a  main-e*fect  plan.  Hov.  Cv^er  for  main-effect  plans 
the  condition  of  equal  frequencies  is  not  a  necessary  one.  We  will  show 
that  a  necessary  and  sufficient  condition  that  the  estimates  of  the  main 
effects  of  any  two  factors  in  a  main-effect  plan  be  uncorrelated  is  rhat  the 
levels  of  one  factor  occur  v4th  each  of  the  levels  of  the  other  factor  with 
proportional  frequencies.  The  condition  of  proportional  frequencies,  will 
be  deduced  for  a  main-effect  plan  on  two  factors,  A  and  B,  occurring  at 
r  and  s  levels,  respectively..  This  waS  stated  first,  it  is  believed,  by 
Plackett  (1946)  but  his  proof  was  found  to  be  obscure.  Therefore  a 
related  proof  is  presented  below. 

If  the  plan  is  orthogonal  then  the  estimate  of  any  component  of  factor 
A  is  orthogonal  with  the  estimate  of  any  component  of  factor  B.  Let  the 
components  of  factor  A  be  represented  by  (r-1)  orthogonal  contrasts, 
and  the  components  of  factor  3  by  (s-1)  orthogonal  contrasts.  Denote 
by  A^  and  B^  the  u-th  orthogonal  contrast  among  the  r  levels  of 
factor  A  and  the  v-th  orthogonal  contrast  among  the  s  levels  cf 
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factor  B,  respectively.  Denote  by  the  coefficients 

of  A  ,  and  by  b,  ,  b..  ,  b  the  coefficients  of  B  .  The  model 

u  ’  Iv  2v  ’  sv  V 

■which  exhibits  these  orthogonal  contrasts  is 


r-1 

s- 

1 

y..  =  a  + 

s 

a.  A  +  S 

b.  B 

ij 

u-l 

lU  u 

1  J'' 

j  =  0.  1, 

2, 

•  •  •  »  (s*  a)  , 

where 

ij 

combination  for  which  factor  A  occurs  at  the  i  level  ara  factor  B 
occurs  at  the  j  level,  p,  is  the  overall  mean  and  e.j  is  the  experi¬ 


mental  error  associated  with  the  observed  yield  y. 


Let  n  =  the  number  of  trials  in  the  plan, 


n.  =  the  number  of  times  the  i  level  of  factor  A  occurs 
in  tue  plan, 

n  =  the  number  of  times  the  j  level  of  factor  B  occurs 
‘  in  the  plan. 


n. .  =  the  number  of  times  the  i  level  of  factor  A  occu.s 
with  the  1  level  of  factor  B. 


Hence 


n. 

1. 


n  . 

•J 


and 


n. 


Theorem  2:  A  necessary  and  sufficient  condition  that  the  estimates  of  the 

components  of  two  factors  A  and  B,  in  a  main-effect  plan,  be  orthogonal 

to  each  other  and  also  to  the  mean  u  is  that  n. .  =  n.  n  ,/n. 

ij  1.  .  J 

Proof:  In  order  that  the  estimates  of  the  components  of  factors  A  and  B 
e  orthogonal  to  each  other  and  also  to  the  nean,  the  design  matrix  X 
must  be  such  tha'  X'X  is  a  diagonal  matrix.  With  the  model 
r-1  s-1 

y..=u+Sa.  A+2b.  B+  e..;  i  =  0,  1,  2,  ....  (r-1) ; 

•’  u=  1  V- 1  ■  •’ 

j  =  0,  I,  2,  ...  ,  (s-1),  the  following  equations  must  hold  in  order  that 
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the  design  matrix  be  of  the  required  form; 


and 


S  a.  n.  =  0  ;  u  *-  1 ,  2,  ♦ . .  , 

lU  l-  *  ’ 

1 

(r-1) 

(5) 

Zb.  n  =0;v=I,  2, 
j  JV  .J 

(s-1) 

(6) 

2  a.  a.  n.  =  0  u  ^  u* 

lU  lU’  l. 
i 

(7) 

2  b.  b.  n  .  =  0  ;  V  ^  v’ 

j 

(8) 

2  a.  b.  n..=  0;u=i,  2,  .... 

.  .  lU  JV  11 

i.J  ^ 

(r-l);v.l.  2.  ....  (s-1). 

(9) 

Equations  (5),  (6)  and  (9)  can  be  expressed  in  matrix  notation  by 
equations  (10),  (11)  and  (12); 


A»N 

r. 


=  d 


■1.1 


(10) 


where  A’  is  an  (r-1)  x  r  matrix  of  coefficients  of  A^ , 

NJ,  ~  (*^0  *^1  ■*’  ^(r-l)  ®mn  ***  nixn  matrix  of  zeros; 


B»N 


.  8 


S-1,  1 


(11) 


where  B*  is  an  (s-1)  x  s  matrix  of  coefficients  of  B^,  and 
^‘s=  Ko^l  •••  ".(s-1)  ’ 


A»NB  =  0 

r 


1,  s- 1 


(12) 


where  N  =  (n.. ). 

'  ir 


Equations  (7)  and  (8)  art  aatumatically  satisfied  since  the  a^^  and  the 

bj^  are  coefficients  of  the  orthogonal  contrasts.  Thus  we  need  only  show 

that  a  necessary  and  sufficient  condition  that  A’NB  =0  ,  ,  .  ei'/en 

r-l, s-1  ^ 
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lhai  A‘i\  =0  ,  ,  and  B’N  =0  ,  , ,  is  that  n..  =  n.  r  ./n, 

r,  r-1,1  .s  ij  1.  .j 

which  expressed  in  matrix  notation  is  N  =  N’^/n. 

To  show  that  this  condition  is  sufficient,  assume  that  N - 


Tren  A’NB  =  A>N  N*  E/u  =  0  .  since  A'N  =Q  ,  , 

r.  .  s  r-i,  s-1  r.  r-1, 1 


tnd 


B  =  »,  , . 

s.  l,s-l 


To  show  that  this  condition  is  also  necessary,  assume  thiit 

A’NB  =  0  ,  ,  -  Since  n.  =  2  n. .  and  n  .  =  2  n. then  N  =  NE  , 

r-I,s-l  1.  ^  ij  .j  ^  ij’  r,  si 

and  M’  -  E,  N,  where  is  an  mxn  matrix  whose  elements  are 

.  s  Ir  mn 

all  \mity.  Let  P  =  :  -^J  ■  Since  the  columns  of  A  are  the 

coefficients  of  (r-l)  orthogonal  contrasts,  P  must  be  non-singular. 
Let  Q  =  I  E^j  ;  bJ  .  Since  the  columns  of  B  are  the  coefficients  of 
(s-1)  orthogonal  contrasts,  Q  must  be  non-singular. 


Now 

P*NQ  = 

1 

.  .  :  B 
si : 

n  N»  b"^ 

.s 

n 

®1,  S-1 

A«N  A*NB 

L.  r-  J 

®r-l,  1 

®r-l,s-l 

Thu*!  P*NQ  is  of  rank  one.  Since  P  and  Q  are  both  non-singular 
matrices,  M  must  have  a  rank  of  one.;  Hence  each  row  of  N  is  a 
multiple  of  the  first  row  and  each  column  is  a  multiple  of  the  first  column. 
Therefore  n../n.  =  n  ./n  or  n. .  =  n.  n  ./n  which  can  be  expressed  in 

ij  i.  .j  ij  1.  .j 

matrix  not.ution  as  N  =  N  N*  /n. 

r.  .  s 

The  theorem  can  easily  be  generalized  to  prove  that  a  necessary  and 
sufficient  condition  that  the  esti.nates  of  the  components  of  k  factors  in 
a  main -effect  plan  be  pairwise  orthogonal  and  also  orthogonal  to  the  mean 
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|i  is  that  the  levels  of  each  factor  occur  with  the  levels  of  any  other 
factor  with  proportional  frequencies.  This  generalization  can  be  made  by 
showing  that  for  any  pair  of  factors  the  proportional  frequency  condition 
is  both  necessary  and  sufficient  to  yield  orthogonal  estimates. 


D.  Plans  for  Asymmetrical  Factorial  Experiments 

Tf  the  levels  of  each  factor  are  arrang-cd  so  that  they  occur  with  the 
levels  of  any  other  factor  with  proportional  frequencies,  it  is  possible  to 
derive  new  classes  of  orthogonal  main-effect  plans  for  asymmetrical 
factorial  experiments.  One  such  class  permits  the  estimation  of  all 
main  effects  without  correlation  for  an  experinic.nt  involving  t^  factors 
at  3j  levels,  t^  factors  ai  s^  levels,  up  to  tj^  factors  at  Sj^  levels, 
with  s^  trials,  where  Sj  is  a  prime  or  the  power  of  a  prime, 


s ,  >  s  -  >  . >  s,  and 
12  k 


2  t  <  (s’?  -l)/(s  -  1). 
i=  I 


A  method  of  constructing  an  orthogonal  main-effect  plan  for  the 
t,  t,  t, 

j  ^  n 

Sj  X  X  ...  X  Sj^  experiment  in  trials  involves  collapsing 
factors  occurring  at  levels  to  factors  occurring  at  levels 
(i  =  2,  3,  4,  ...»  k)  by  utilizing  a  many-one  correspondence  of  the  set 
of  Sj  levels  to  the  set  of  s.  levels.  First  construct  an  orthogonal 
main-effect  plan  for  the  symmetrical  factorial  experiment  involving 
(Sj  -i)/(Sj-l)  facto  s,  each  at  levels,  with  s”  trials,  where  Sj 
is  a  prime  or  the  power  of  a  prime.  Collapse  the  levels  of  t^  of  these 
faetD's  to  s,,  levels,  where  s.,  <  s.,  by  making  a  many-one 

it  1 
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correspondence  of  the  set  of  levels  to  the  set  of  s^  levels. 

Similarly  collapse  the  levels  of  t^  of  the  original  factors  to  levels, 

where  s- <  s^  <  s,,  and  so  on. 

If  for  some  i,  s,  -  ,  then  a  factor  with  s,  levels  can  be 

collapsed  into  (s^-l)/(Sj^-l)  factors  each  having  Sj,  levels.  Since  there 

exists  an  orthogonal  main-effect  plan  for  (s^ -l)/(Sj^-l)  factors,  each 

at  s.  levels.,  with  s^  treatment  combinations,  we  can  replace  each  of 
11 

the  Sj  levels  by  one  of  the  s^  =  s^  treatment  combinations.  To 
illustrate  this  pcini  consider  a  factor  at  s^  =  4  levels.  There  exists  an 
orthogonal  main-effect  plan  for  three  factors,  each  having  two  levels,  in 
four  treatment  combinations,  namely:  00  0,  Oil,  101  and  110. 

If  we  make  the  following  correspondence: 

Fo'-ir-level  Two-level 

factor  factors 

0  - »  0  0  0 

1  - !•  0  1  1 

2  - *  1  0  1 

3  »  1  1  0 

the  four-level  factor  is  collapsed  to  three  two-level  factors. 

If  the  (s.-l)  degrees  of  freedom  for  each  of  the  t.  factors  at  s. 
levels  are  represented  by  (s.-l)  orthogonal  contrasts  among  the  s. 
levels,  the  estimates  of  these  contrasts  for  any  factor  will  be  uncorre- 
latcd  with  the  estimates  of  the  contrasts  for  any  other  factor  because  the 
Correspondence  scheme  automatically  guarantees  proportional  frequencies 
of  the  levels  of  each  factor. 

2  2 

An  orthogonal  .main- effect  plan  for  the  2x3  factorial  experiment 
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with  nine  trials  is  now  constructed  to  illustrate  the  technique  of 


collapsing  levels.  First  construct  an  orthogonal  main-effect  plan  fo 


four  factors,  each  having  three  levels  with  nine  treatment  combinations. 

0  0  0  0 

Oil?. 

0  2  2  1 
10  11 
112  0 
12  0  2 
2  0  2  2 
2  10  1 
2  2  10 

Collapse  each  of  the  first  two  factors  to  two-level  factors  using  the 


followinp  correspondence  scheme; 


Three-level 

factor 


0 

1 

2 


Two-level 

factor 

0 

1 

0 


The  resulting  treatment  combinations  constitute  an  orthogonal  main- 


2  2 

effect  plan  for  the  2  x  3  experiment  and  are  displayed  below. 


0 

0 

0 

0 

1 

0 

0 

2 

f\ 

1 

I 

* 

2 

0 

0 

2 

2 

0 

c 

2 

1 

0 

1 

0 

1 

1 

c 

1 

1 

0 

0 

1 

0 

1 

1 

2 

0 

Doubling  the  number  of  trials  and  doubling  the  number  of  levels  of  one 


factor  leads  also  to  some  new  orthogonal  main-effect  plans.  To 

4 

illustrate  the  construction  procedure  consider  the  3  plan  in  9 


observations  and  repeat  it,  replacing  the  levels  0,  1  and  2  in  one  of 
the  factors  by  the  levels  3,  4  and  5.  This  gives  a  6  x  3^  plan  ’.n  18 
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'*1  ”2 

trials.  The  collapsing  procedure  will  then  give  the  5x3  x  2  plan, 


S  n.  -  3,  in  18  trials, 
i.l  ^ 


1  ne  Class  oi  ortnot'c 


f  .1.  1  2  -k 

ins  * —  3  X...XS, 

i  2  k 


experiment  with  s,  trials  where  s,  >St  >..._>  s,  irestricts  the  number 
of  trials  to  be  equal  to  s^  where  is  the  lajrgi-est  number  of  levels.. 
Thus,  for  example,  one  would  require  sixteen  trials  in  order  to  construct 

4 

an  orthogonal  main-effect  plan  for  the  4x2  experiment  using  the 
procedures  suggested  above.  A  second-class  of  orthogonal  main-effect 


’  I  “'2  K 

plans  can  be  derived  for  the  s.  xs,  x...  xs  experiment  in  s 


1 


trials,  where  Sj^  is  a  prime  or  the  power  of  a  prime,  s^^tt  Sj^, 

k 

and  Z  X.  t.  <  (s*^  -  l)/(3,  -  1)  where  1  =  X.<X,<A,<  ...<X,,  the  X. 

1=  1 

being  integers.  An  orthogonal  main-effect  plan  of  this  class  exists  for 

4 

the  4x2  experiment  with  only  eight  trials. 


Theorem  3:’  Consider  an  orthogonal  main-eflEect  plan  for  (s*^  -  l)/(s  -  i) 

factors,  each  at  s  levels,  with  .s^  trials,  where  s  is  a  prime  or  the 

power  of  a  prime  number.  Then  a  factor  at  t  levels,  where 
2 

s  <t  <  s  ,  can  be  introduced  as  a  replacement  .for  a  suitably  chosen  set 
of  (s  ¥  1)  factors  in  such  a  \  ly  as  to  preserve  orthogonality  of  main- 
I'ffect  estimates. 


Proof:  Let  t  =  .  There  exists  an  ortnogonal  main-effect  plan  for 

2  2 
(&  -  l)/(s  -  I)  factors  each  at  s  levels  in  t  =  s  trials.  Hence  a  factor 

2  2 

having  t  =  s  levels  can  replace  (s  -  l)/(s-  1)  =  (a  +  I)  factors  each 


30 


having  s  levels.  D  t  <  s  then  a.  factor  having  s”  levels  can  replace 
(5+  1)  factors  each  having  s  levels  and  then  collapsed  to  a  t-ievel 
factor  by  a  many-one  correspondence  scheme. 

2 

Corollary  1:  The  maximum  nambex  of  t-level  factors  (s  <  ■■- <  s  ) 
which  can  be  introduced  into  an  orth(^cnal  main^.«£ect  plan  fer 
{s^-l)/(s-l)  factors,  aach  at  s  levet^.  with  trials  is 
(i)  {s^-l)/(s^-l)  if  n  is  even  and  (ii)  the  largest  integer  less  than  or 
equal  to  f  -  IJ  tf  n  is  odd. 

1  rH 

Corollary  2:  A  factor  at  t  levels,,  where  s'  <  t<  d  can  be 

introduced  as  a  replaceresent  for  a  suitably  ch  issen  set  of  l)/(s-l) 

factors  each  having  s  levels  in  such  a  way  as  to  jKreserye  the 
orthogonality  of  main- effect  estimate®. 

This  replacement  procedure  wtll  be  iHust^-ated  by  construct  -ig  an 

A 

orthogonal  main- effect  plan  for  the  4x2*  experiment  in  eight  trial.',. 
First  construct  '"1  orthogonal  naain-effecC  plan  the  2*  ex{>erimeBt. 
The  seven  two-level  factors  caa.i>e  represented  by  SEj.,.  X^.  Xj  +  X^, 
X^,  Xj  T  X,,  X^  i  X^  and  X^  +  X^  +  X^  ..  The  treatment  combinations 
for  this  plan  are  the  fcBtowingt 

0  0  0  0  0  e  0 

0  0  0  till 
0  I  1  0  0  1  1 
0  1  1  1  1  0  0 
1  0  1  0  I  0  I 
1  0  I  I  0  1  0 
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110  0  110 


110  10  0  1 

It  is  knovm  that  there  exists  an  orthogonal  main-effect  plan  for  the 
exj^rimenr  sa  four  trials.  The  treatment  combinations  for  this  pla" 
a-j?e  yt  #  uj.,  (0  1  1),  fl  0  1)  and  (1  t  0).  Thus  by  choosing  three 
faenors  a£  the  2^  plam  whose  X  representations  are  such  that  the 
genes^ifficii  iate^iaction  of  any  two  of  the  three  factors  is  the  third 
factory,  two— ievel  fee*  can  be  replaced  by  a  four-level  factor 

accordag-tis-iftwr  following  correspondence  scheme 


T  wo—iev-el 

F  our -level 

fectors 

factor 

0  C  0 

- > 

0 

0  1  1 

- ^ 

1 

1  0  1 

- :> 

2 

1  I  0 

- 

3 

Since  the  X  ’•epresentations  of  fiie  Sxst  three  factors  of  the  above  plan 
are  X,,  and  Xj  +  X^,  these  three  factors  can  be  replaced  by  e 

4 

;our-level  factor  and  the  orthogonal  main'flffect  plan  for  the  4x2 
experiment  in  eight  trials  is  given  by  the  following  treatment  combinations 

0  0  0  0  0 

0  till 

1  0  0  11 

1  110  0 

2  0  10  1 

2  I  ®  I  0 

S2 


3  0  110 


3  10  0  1 

By  collapsing  the  four-level  factor  to  a  three-level  factor,  an 

4 

orthogonal  main-effect  plan  for  the  3  x  Z  experiment  is  obtained. 

It  can  be  easily  verified  that  a  suitably  chosen  set  of  {s^-l)/(s-f) 

3 

factors,  each  having  s  levels,  occurring  in  an  orthogonal  main-effect 
6  3 

plan  with  s  trials  can  Le  replaced  by  (s  -l)/{s-l)  factors,  each 
2 

having  s  levels.  This  p-  jposition  can  be  illustrated  by  replacing 
three  eight-level  factors  by  seven  four-level  factors  in  an  orthogonal 
main- effect  plan  with  sixty-four  trials. 

63 

Consider  the  orthogonal  main-effect  plan  for  the  2  experiment  in 
sixty-four  trials.  Let  each  factor  be  represented  by  an  effect  or 
interaction  of  the  2^  factorial  experiment,  namely  Xj,  X^,,  X^,  X^, 

X.,  X^  or  any  one  of  their  generalized  interactions.  From  Corollary  2 
of  Theorem  3  it  is  clear  that  each  eight-level  factor  introduced  into  the 
plan  replaces  seven  two-level  factors.  Let  us  denote  three  eight-level 
factors  by  A,  B  and  C.  Table  12  gives  the  X  representations  for  the 
two-ievel  factors  which  -are  replaced  by  the  three  eight-level  factors. 

It  will  be  noted  in  Table  12  that  the  X  representations  of  the  two- 
level  factors  which  are  replaced  by  the  eight-level  factor  C  are  the 
general’ red  interactions  of  the  X  representations  of  the  two-level 
factors  which  are  replaced  by  factors  A  and  B.  Thus  each  row  of  the 
table  represtnts  three  two-level  factors  which  can  be  replaced  by  a 
four-level  ..actor..  Hence  it  is  clear  that  three  eight-level  factors  can  be 
replace!"  by  seven  four-level  factors. 

G3 


TABLE  12 


TWO-LEVEL  FACTORS  REPLACED  BY  EIGHT-LEVEL  FACTORS 


A 

B 

C 

^1 

X1+X2 

X^rXs 

Xj  +  X^  +  X^+Xg 

X1  +  X2+X5 

x^  +  x, 

4  -3 

X3-X3 

Xj-.  x^+x^  +  x.+x^ 

X^4X,4X, 

XJ+X3+X5 

X,+X34X5  +  X^ 

^6 

x.,  +  x,  +  x*  +  x. 

2  3  4  5 

X2^X34X^ 

Xj  tX^+Xj  +  X^  +  Xg 

^3  +  ^6 

Xi  +  X^  +  X^  +  X.  +  X^ 

It  is  evident  that  the  use  of  factors  for  which  the  levels  occur  with 
proportional  frequencies  also  yields  orthogonal  main-effect  plans  for 
symmetrical  factorial  experiments.  For  example,  an  orthogonal  main- 

5 

effect  plan  for  the  3  experiment  can  be  constructed  with  sixteen  trials 

5 

by  collapsing  all  the  four-level  factors  in  the  plan  for  the  4  experiment 
to  three-level  factors. 

The  use  of  factors  whose  levels  occur  with  proportional  frequencies 
also  peiTnits  the  construction  of  orthogonal  main-effect  plans  for  factors 
for  which  the  number  of  levels  is  not  equal  to  a  prime  or  the  power  of  a 
prime.  One  such  plan  with  forty-nine  trials  permits  uncorrelated  main- 

g 

effect  estimates  for  the  6  experiment.  This  plan  can  be  const r  Acted  by 
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''ollapsing  the  seven-level  factors  to  six-level  factors  in  the  plan  for  the 


7  experiment. 


E.  Efficiencies  of  Main-Effect  Estimates 


Although  any  many-one  correspondence  of  the  set  of  s^  levels  to  the 
Set  of  s^  levels  will  yield  proportional  frequencies  of  the  levels,  there 
arises  the  piobJem  of  which  correspondence  is  "best"  in  some  sense. 

The  problem  may  be  solv-  a  by  determining  the  efficiencies  of  the  main- 
effect  estimates  obtained  usitig  proportional  frequencies  relative  to  the 
estimates  which  would  result  rom  using  equal  frequencies  of  the  levels 
of  each  factor. 

As  an  illustration  we  will  calculate  the  relative  efficiency  of  a  three- 
level  facte  r  in  a  main-effect  plan  with  twenty-five  trials. 

Assume  the  correspondence  scheme  used  to  collapse  a  five-level 
factor  to  three  levels  is  as  follows; 


Five-level 

factor 


Three-level 

factor 


The  levels  0,  1,  and  2  occur  in  the  ratio’s  2:1:2.  Thus  for  this 
factor  the  0  level  occurs  in  ten  treatment  combinations,  the  1  level 
occurs  in  five  treatn*eiit  combinations  and  the  2  level  occurs  in  ten 


treatment  combinations. 
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The  variance  of  the  linear  effect  estimate  of  this  factor  is  equal  to 

2 

cr  /20  and  hence  the  information  on  a  unit  basis  is  equal  to 

20/25 O'"  =  4/5  (T^.  The  variance  of  the  linear  effect  estimate  of  a  three- 

Icvcl  factor  in  3  trials  is  equal  to  ^  /2*  3  and  the  inTormatlon  Oh  a 

unit  basis  is  2.  3*^  =  2/3<r^.  Hence  the  relative  efficiency  of  the 

linear  effect  estimate  is  equal  to  4/5  x  3/2  =  6/5. 

The  variance  of  the  quadratic  effect  estimate  for  the  three-level 

2 

factor  in  twenty-five  trials  .s  equal  to  <r  /4  and  the  information  is  then 

equal  to  4/25  0-^.  The  variance  of  the  quadratic  effect  estimate  with  3*^ 

trials  is  equal  to  (r^/2.  3*^  ^  and  hence  the  information  on  a  unit  basis  is 
2 

equal  to  2/9  <r  .  The  relative  efficiency  of  the  quadratic  effect  estimate 
is  therefore  equal  to  4/25  x  9/2  =  18/25. 

The  relative  efficiencies  of  the  estimated  effects  are  presented  for 
various  proportional  frequencies  in  Table  13  .  One  would  chose  the 
proportional  frequencies  which  give  the  greatest  efficiency  of  estimates. 
Thus  for  example,  if  an  experiment  in  twenty-five  trials  involved  two-levcl 
factors  the  two  levels  should  occur  in  the  ratio  2:3  rather  than  in  the 
ratio  1  :  4  because  the  efficiency  of  the  2  :  3  ratio  is  24/25  whereas 
the  efficiency  of  the  1  :  4  ratio  is  only  16/25. 
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TABLE  13 


RELATIVE  EFFICIENCIES  OF  MAIN-EFFECT  ESTIMATES 


Level 

0  1 

Efficiency 

Proportional  frequency 

1  :  2 

8/9 

Z  :  3 

24/25 

1  :  4 

16/25 

3  ;  4 

48/49 

2  :  5 

40/49 

I  :  6 

24/49 

Level 

0  I  2 

Contrast 

Proportional  frequency 

Linear 

1  ;  2  :  i 

3/4 

Quadratic 

1:2:1 

9/8 

Linear 

2:1:2 

6/5 

Quadratic 

2:1:2 

:.8/25 

Linear 

1:3:1 

3/5 

Quadratic 

1:3:1 

27/25 

Linear 

2:3:2 

6/7 

Quadratic 

2:3:2 

54/49 

Linear 

3:1:3 

9/7 

Quadratic 

3:1:3 

27/49 

Linear 

1:5:1 

3/7 

Quadratic 

1:5:1 

45/49 
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F.  Blocking  in  Main-Effect  Plans 


Even  though  the  orthogonal  main-effect  plans  art  highly  fractionated 
these  plans  may  still  require  more  trials  than  can  he  carried  out  under 
uniform  conditions.  Thus  it  would  be  desirable  to  divide  the  experi¬ 
mental  data  into  smaller  blocks  in  such  a  manner  that  the  main  effects 
may  still  be  estimated  without  correlation.  In  this  section  we  will 
illustrate  the  use  of  some  of  the  factors  in  an  orthogonal  main-effect 
plans  as  blocking  factors. 

4 

Consider  the  orthogonal  main-effect  plan  for  the  3  experiment  with 
nine  trials.  The  treatment  combinations  for  this  plan  are 

0  0  0  0 
0  112 
0  2  2  1 
10  11 
112  0 
12  0  2 
2  0  2  2 
2  10  1 
2  2  10 

If  there  are  only  three  factors  under  investigation  the  fourth  factor  of  the 
above  plan  can  be  used  as  a  blocking  factor  to  yield  the  following  blocks: 


Block  1 

Block  2 

Block  3 

0  0  0 

0  2  2 

0  1  1 

1  1  2 

1  0  1 

12  0 

2  2  1 

2  1  0 

2  0  2 

S8 


The  esiimate  of  the  main  effects  of  the  three  factors  are  clear  of  the 
biock  effects  since  each  level  of  the  three  factors  occurs  once  in  each 
block.  The  linear  effect  of  the  first  factor  is  given  by 

g  {2  21  +  210+202  -  0  00-022-01  1). 

It  is  evident  that  each  block  effect  enters  this  estimate  once  positively 
and  once  negatively  and  hence  the  estimate  is  clear  of  block  effects. 

If  the  four  factors  in  the  orthogonal  main-effect  plan  are  represented 
by  Xj,  X^,  Xj  +  X2  and  +  2X2  the  use  of  the  fovirth  factor  as  a 
blocking  factor  is  equivalent  to  confounding  the  factor  represented  by 
X j  +  2  X2  with  blocks. 

In  general,  if  two  factors  are  used  as  blockine  factors  then  so  are  the 

factors  represented  by  the  generalized  interactions  of  their  X  represen- 

7 

tations.  For  example,  if  the  seven  factors  in  the  plan  for  the  2 
experiment  with  eight  trials  are  represented  by  Xj,  X21  Xj  +  X2> 

Xj,  Xj  +  Xj,  X2+X.^  and  XJ  +  X2+X2  an  orthogonal  main-effect  plan 

A 

for  the  2  experiment  in  4  blocks  of  2  treatment  combinations  can  be 
obtained  by  using  the  factors  represented  by  Xj,  X21  Xj  +  X2  as 
blocking  factors.  This  is  equivalent  to  confounding  Xj,  X2  and  X1  +  X2 
with  blocks. 

2  6 

Now  consider-  the  orthogonal  main-effect  plan  for  the  4x3x2 
experiment  with  yixteen  trials.  The  plan  is  comprised  of  the  following 
treatment  combinations: 
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000  000000 

oil  lOiilO 

022  110011 

fill  0  1110  1 

10  1  0  110  11 

110  110  10  1 

121  101000 

112  000110 

202  101101 

211  000011 

2  2  0  0  1  1  1  1  0 

2  11  1  1  0  0  U  0 

301  110110 

3  1  2  0  1  1  0  0  0 

321  000101 

3  10  10  10  11 


The  following  orthogonal  main-effect  plans  which  utilize  one  or  more 
of  the  factors  as  blocking  factors  may  be  constructed  from  the  given  plan: 
3  5  • 

(i)  4x3x2  in  2  blocks  of  8  treatment  combinations; 

Using  the  last  two-level  factor  as  a  blocking  factor  the  two  blocks 
consist  of  the  treatment  combinations  presented  below: 
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Block  I 


Block  Z 


0  0  0 

0  0  0  0  0 

0  2  2 

110  0  1 

0  1  1 

10  111 

0  1  1 

0  1110 

1  2  1 

10  10  0 

1  0  1 

0  1  1  0  1 

112 

0  0  0  1  1 

1  1  0 

110  10 

2  2  0 

0  1111 

2  0  2 

10  110 

2  1  1 

110  0  0 

2  1  1 

0  0  0  0  1 

3  0  1 

110  11 

3  2  1 

0  0  0  1  0 

3  1  2 

0  110  0 

3  10 

10  10  1 

2  3 

(ii)  4x3  X  2  in  4  blocks  of  4  treatment  combinations: 

Consider  the  last  three  two-level  factors  only.  The  levels  for 
these  three  factors  occur  in  the  four  sets  00  0,  Oil,  101  and 
110  each  occurring  four  times  in  the  sixteen  trials.  Th  treat¬ 
ment  combinations  of  the  first  six  factors  can  be  blocked  according 
to  the  paiticular  set  to  which  the  levels  of  the  last  three  factors 
belong.  Hence  the  four  blocks  are: 


Block  1 


Block  2  Block  3 


Block  4 


{iii> 


0  0  0  0  0  0 
12  110  1 
2  11110 
312011 


022110 
10  10  11 
2  110  0  0 
3  10  10  1 


0  110  11 
110  110 
2  0  2  1  0  1 
3  2  1  0  0  0 


X  2  in  4  blocks  of  4  treatment  combinations: 


0  1110  1 
1  1  2  0  0  0 
2  2  0  0  1  1 
3  0  1110 


Utilizing  the  four-level  factor  as  the  blocking  factor  the  four 


blocks  are: 
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Block  I 
OCOGOOOO 
i  1  1  0  1  1  I  0 
2,  2  1  1  0  0  1  1 
110  1110  1 
,6  . 


Block  2 
0  10  110  11 
10  110  10  1 
21101000 
12000110 


Block  3 
0  2  10  110  1 
11000011 
20011110 
11110000 


Block  4 
0  1110  110 
12011000 
2  1  0  0  0  1  0  1 
J.  0  I  0  1  0  1  1 


(iv)  4x3x2  in  4  blocks  of  4  treatment  combinations: 

This  plan  can  be  constructed  by  considering  the  first  three-level 

factor  tc  be  a  foui  level  factor  and  using  that  factor  as  a  blocking 

factor.  If  every  second  1  in  the  first  three-lcvcl  factor  of  the 

2  ^ 

main-effect  plan  fcr  the  4x3x2  experiment  is  replaced  by  a 
3  the  sixteen  treatment  combinations  then  comprise  a  main-effect 
plan  for  the  4^  x  3  x  2^  o^'eriment. 


Block  1 


Block  2  Block  3 


OOOOOOOO 
110  110  11 
ZZlGiiOl 
3  1110  110 


0  110  1110 
10  110  10  1 
21000011 
32011000 


0  2  1  1  0  0  1  1 
11101000 
20011110 
31000101 


Block  4 
0  10  1110  1 
12000110 
21110000 
30101011 


It  is  clear  that  in  each  of  the  above  plans,  the  main-effect  estimates 
and  the  block  effect  estimates  are  uncorrclated. 
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G.  Randomizaticu  Procedure 


An  important  aspect  of  most  experimental  situations  is  the  fact  that 
each  experimental  unit  can  be  subjected  to  only  one  of  the  treatments  of 
interest.  Because  of  this  fact  the  variability  due  to  heterogeneity  of 
experimental  units  w.li  contribute  to  experimental  uncertainty.  To  obtain 
some  control  of  this  variability  the  device  of  randomization  is  used  in 
the  statistical  design  of  experiments.  This  technique  implies,  essentially, 
that  random  methods  of  selection  and  assignntent  are  employed  in 
carrying  out  the  experiment. 

The  procedure  recommended  for  assigning  treatments  at  random  to 
the  experimental  units  of  an  orthogonal  main-effect  plan  is  as  follows: 

(i)  Choose  the  appropriate  plan. 

(ii)  Randomly  assign  the  factors  to  the  columns  of  the  chosen  plan. 

(iii)  Randomly  assign  the  levels  of  each  factor  to  the  numbers 
0,  1,  2,  ...»  representir.g  the  levels  of  a  factor. 

(iv)  Randomly  assign  the  treatments  to  the  experimental  unit?. 

To  illustrate  this  procedure  we  shall  describe  the  randomization 
procedure  to  be  followed  with  an  experiment  involving  three  factors 
A,  B  and  C.  each  having  three  levels  and  one  factor,  D,  at  two  levels. 
The  appropriate  orthogonal  main-effect  plan  for  this  experiment  is  given 
by  the  following  nine  treatment  combinations. 
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q  o 


VH  1  [,  A  \;w  A  ±'  3  A  V  .  S. 


0  0  0  0 

0  110 

0  2  2  1 

10  11 
112  0 

1  2  G  0 

2  0  2  0 

2  10  1 

2  2  10 

Assign  the  factors  A,  B  and  C  at  random  to  the  first  three  columns 
of  th-  above  plan  and  assign  factor  D  to  the  fourth  column.  Then,  for 
each  of  the  factors  A,  B  and  C  randomly  assign  the  three  levels  to  0, 

1  and  2.  Similarly  for  factor  D  assign  the  two  levels  to  0  and  1  at 
random.  Then  these  treatments  are  assigned  to  nine  experimental  units 
at  random,  for  example,  by  testing  the  combinations  in  random  order. 

H.  Analysis  of  Main-Effect  Experiments 

At\  important  feature  of  the  full  factorial  arrangement  is  that  the 
main  effects  and  all  interactions  can  be  estimated  without  correlation 
Since  the  main-effect  plans  developed  in  this  report  allow  uncorrelated 
estimates  of  all  main  effects  the  analyses  of  these  experiments  are 
similar  to  the  analysis  of  a  complete  factorial  experiment.  Estimation 
is  based  on  the  general  procedure  described  in  Chapter  II,  and  a  quick 
review  of  aspects  relevant  to  main  effect  plans  will  now  be  given. 

The  multiple  regressio.:  model  for  an  orthogonal  main-effect  exp*‘ri~ 
ment  can  be  written  in  matrix  notation  as  y  =  Xp  +  e  where  p  is  the 
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vector  of  effects  and  interactions.  The  estimates  of  the  effects  and 
interactions  are  given  by  p  =  (X*X)  ^  X’y,  v/liere  (X*X)  ^  is  the 
variance-covariance  matrix.  The  property  of  uncorrelated  estimates  is 
reflected  in  the  fact  that  the  variance-covariance  matrix  is  a  diagonal 
matrix. 

To  illustrate  the  estimation  procedure  we  consider  the  plan  for  two 
two-level  factors,  A  and  B,  and  two  three-level  factors,  C  and  D, 
the  levels  being  equally  spaced,  in  nine  trials,  when  all  interactions  are 
assumed  to  be  absent.  The  plan  is  giver,  by  the  following  set  of  treatment 
combinations: 

0  0  0  0 

0  112 

0  0  2  1 

10  11 
112  0 
10  0  2 
0  0  2  2 

0  10  1 

0  0  10 


The  responses  may  be  expressed  in  terms  of  the  main  effects 


as 


y...  -  u  +  3..  A  b.  B  +  c.  G,  e  C,..,  +  d  D-  +  d  D.  +  e.., 

■'ijkm  '^1  j  li  K,,  Q  m^  L  m,..,  Q  ijkm 


Q 


Q 


where  A,  B,  Cj^,  Cq,  and  Dq  are  the  effects  of  the  respective 

factors  and  a.,  b. ,  c,  ,  c,  ,  d  ,  d 

i’  1  ’  k  '  k_  ’  m.  m^ 

•’  L.  Q  L  Q 


are  the  coefficients  of  the 


orthogonal  contrasts  defining  the  corresponding  effects.  The  factors  are 
assumed  to  be  quantitative  factors  and  the  levels  of  the  factors  are 
equally  spaced.  * 

If  any  of  the  factors  are  qualitative,  they  can  still  be  treated  as 
quantitative  factors,  except  that  what  are  contrasts  of  specific  meaning 
in  the  quantitative  case,  such  as  linear  and  quadratic  effects,  are  merely 
contrasts  among  the  ’evels  of  the  qualitative  factors.  For  example,  if 
we  use  the  numbers  0,  1  and  2  to  denote  the  levels  of  a  factor,  F,  at 
three  levels  and  get  what  look  superficially  to  be  linear  and  quadratic 
effects,  they  are  in  fact 

L  =  F,  -  Fo 

Q  =  F2  -  2Fj  -  Fq  =  (F2  -  Fj)  -  (Fj  -  Fq) 

where  L  and  Q  denote  the  linear  and  quadratic  effects  and  F.  denotes 
the  treatment  combinations  which  contain  factor  F  at  the  i  level.  From 
such  calculated  effects  one  can  determine  any  contrasts  which  seem 
relevant.  For  instance 


F2  -  Fj  =  (L  +  Q)/2 

and 

Fj  -Fo=(L-Q)/2. 


If  the  levels  of  quantitative  factors  are  not  at  equally  spaced  intervals 
the  effects  can  still  be  written  in  terms  of  orthogonal  contrasts.  A 
procedure  for  obtaining  orthogonal  polynomials  for  unequally  spaced  levels 
is  given  in  section  C  of  Chapter  V. 
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T 


The  matrix  of  known  coefficients*  is  given  by 


(A 

3'"Q 

I°Q 

1 

-1 

-1 

-1 

1 

-1 

1 

1 

-1 

2 

0 

-2 

1 

1 

1 

-1 

-1 

1 

1 

0 

-2 

X  = 

1 

2 

-1 

0 

-2 

n 

-2 

1 

2 

2 

1 

1 

-1 

1 

1 

2 

-1 

-1 

1 

1 

1 

1 

-1 

-i 

3 

1 

1 

1 

1 

-1 

2 

•t 

-  1 

1 

0 

_2 

LI 

-1 

-1 

0 

_2 

-1 

1 

the  information 

matrix 

is 

• 

9 

0 

0 

0 

0 

0 

0 

0 

18 

0 

0 

0 

0 

0 

0 

0 

18 

0 

0 

0 

0 

X»X  = 

0 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

18 

0 

0 

0 

0 

0 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

18 

The  coefficients  of  the  parameters  are  obtained  by  a  convention 
which  is  discussed  in  section  C  of  Chapter  V. 
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Since  X*X  is  a  diagonal  matrix  the  plan  is  orthogonal  and  the  variance- 
covariance  matrix  is  given  by 


The  estimates  of  effects  and  interactions  ari 


where  y^,  •  •  • ,  y^  are  the  responses  of  the  nine  treatment 

combinations  in  the  order  presented  in  the  plan. 

Thus,  for  example. 


and 


jA  =  A  f  _  y^  +  2 +  Zy ^  +  2y^  -  Yy  -  Vs  "  • 

=  i  ^-yi^ys^ys-ye^yi-ys^  ■ 

The  variances  of  the  estimates  are  obtained  from  the  variance- 
covariance  matrix.  Thus, 


var  (A) 

=  \’ar  (B: 

=  0-^/2 

JK 

=  var(Dj^) 

=  0-^/6 

2 

var  (Cq  ) 

=  varCD^) 

=  a/2 

An  unbiased  estimate  of  a-"  is  derived  from  the  sum  of  squares  of 
deviations  about  the  estimated  values,  namely 

P  =  |(y»y  -P’X^y) 


The  sum  of  squares  in  the  analysis  of  variance  associated  with  any 
contrast  is  merely  the  square  of  the  contrast  divided  by  the  sum  of 
squares  of  the  coefficients  of  the  contrast.  Hence,  the  sum  of  squares 
due  to  is 

r-Yi  -y^  -y3+  2y4+  2y3+  2y^  -y^  -yg  -ygj2 


The  sum  of  squares  due  to  A  is  then 


-  Vz  "  >3  +  -^4  ^''6  ■  y?  -  'S  ■  y' 


.J" 


If  the  total  sum  of  squares  is  corrected  for  the  mean,  the  partitioning 


for  the  analysis  of  variance  is  given  in  Table  14. 


TABLE  14 


PARTITIONING  OF  ANALYSIS  OF  VARIANCE 


Source 

Degrees  of  Freedom 

A 

1 

B 

1 

^L 

1 

1 

Dr 

1 

D 

1 

Error 

2 

Total 

8 

It  is  clear  that  the  two  degrees  of  freedom  available  for  estimation  of 
error  are  the  result  of  collapsing  twc  three-level  factors  to  two-level 
factors.  The  estimate  of  error  can  be  partitioned  into  single  degrees  of 
freedom  as  follows.  Consider  the  levels  of  factor  A  and  factor  B  as 
they  were  before  being  collapsed.  The  levels  are  then  given  as 

Levels  oi  factor  A:  000111222 

Levels  of  factor  B:  0  12012  Oli 

In  order  to  collapse  a  three-level  factor  to  a  two-level  factor  we  make 
the  correspondence 


ICD 


Two-level 

factor 


Three-level 
facUir 

0  - ?  0 

1  - >  1 

2  - >  0 

If  factors  A  and  B  were  three-level  factors  then 

i  ^  ■  ^1  ■  >'2  -  ^3  y?  ^  ^8  ^9^ 

and  =  ^  (  -  y,  +  y3  -  +  Yg) 

Since  these  two  factors  have  only  two  levels  and  the  level  of  factor  A  is 
0  for  each  response  in  the  c  c  A,  and  the  level  of  factor  B  is 

JLj 

0  for  each  response  in  the  estimate  then  these  contrasts  are 

estimating  pure  error.  Thus,  the  two  single  degrees  of  freedom 
Z 

estimates  of  <r  are  given  by 

^  {-yj  -  y2-y3  +  y7  +  y8  +  y9>^ 

and  F^-yi-^y3‘y4‘^y6"y7‘‘y9^^ 

The  partitioning  of  the  analysis  of  variance  is  presented  in  Table  15. 

If  several  estimates  of  error  are  possible  one  can  determine  whether 
they  are  homogeneous  estimates  of  error  (e.  g.  Bartlett’s  test)  and  if 
they  are  found  to  be  hor-.ogeneous  they  can  be  combined  to  give  a  pooled 
estimate  of  error.  Evidence  of  estimates  of  error  being  not  poolable  is 
evidence  that  there  are  interactions  present  in  the  situation,  and  further 
experimentation  to  explore  these  is  needed. 
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TABLE  15 


PARTITIONING  OF  ANALYSIS  OF  VARIANCE 


Source 

Degrees  of  Freedom 

A 

1 

B 

1 

^L 

1 

S 

1 

°L 

1 

1 

Error  A 

1 

Error  B 

1 

Total 

8 
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V.  CATALOGUE  OF  ORTHOGONAL  MAIN-EFFECT  PLANS 


A.  Construction  of  Basic  Plans 

The  task  of  presenting  a  catalogue  which  gives  every  possible 
orthogonal  main-effect  plan  with  81  trials  or  fewer  is  enormous  and  need 
not  he  undertaken.  Each  of  these  plans  can  be  easily  deduced  from  one 
of  twenry-six  "basic  plans",  by  choosing  a  suitable  set  of  columns. 

4 

Consider  the  orthogonal  main-effect  plan  for  the  3  experiment  in 

nine  trials.  It  was  demonstrated  in  Chapter  IV  that  Irom  this  plan  one 

3  2  2,3 

can  obtain  plans  for  the  following  experiments:  3  x2,  3  x2,3x2 

4  4 

and  2  .  If  a  plan  which  consisted  of  the  plans  for  both  the  3  experl- 

4 

ment  and  the  2  experiment  in  nine  trials  is  given,  then  the  plans  for 
^3  2  2*  3  4 

any  one  of  the  3  ,  3  x2,  3  x2,  3x2  or  2  experiments  can  be 

obtained  by  selecting  the  appropriate  number  of  colirmns  from  the  plan 
4  4 

for  the  3  and  2  experiments,  respectively.  The  plan  which  consists 

4  4 

of  the  plans  for  both  tlie  3  experiment  and  the  2  experiment  tn  nine 

trials  is  called  a  basic  plan. 

*3 

Similarly  a  basic  plan  for  the  4  x  3  x  2  e:<pcriment  in  sixteen 

5  5  15 

trials  is  a  plan  consisting  of  the  plans  for  the  4  ,  3  and  2  experi¬ 
ments  . 

Each  column  of  the  basic  plan  represents  a  factor.  The  number  of 
levels  of  any  factor  can  be  determined  by  counting  the  number  of 
different  symbols  0,  1,  2,  etc.  which  represent  the  levels.  The  columns 
are  numbered  so  that  each  column  may  be  identified  quickly.  The 
numbering  of  the  columns  may  best  be  explained  by  an  example.  The 
column  numbers  on  the  four-level  factor®  of  basic  plan  5,  which  consists 
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-5  5  15 

of  the  plans  for  the  4,3  and  2  experiments  in  sixteen  trials,  are; 


1  2  3  4  5. 

¥  *  ♦  * 


The  numbers  on  the  three-level  columns  for  this  plan  are  also 

1  2  3  i  5  . 

=i!  ♦  * 

The  numbers  on  "-he  two-level  factors  range  from  1  to  15  where  for 
tabular  convenience  these  numbers  are  written  in  the  form 

0  0  ....  1  . 

1  2  5 

The  footnote  given  below  the  basic  plan  Indicates  that  the  column 
identified  by  ^  replaces  the  columns  identified  by  ?  ^  and  »  ,  the 

column  identified  by  ^  replaces  the  columns  identified  by  ^  ^  and  ^ 

and  so  on.  Kence,  If  a  four-level  factor  identified  by  ^  Is  used  in  an 

orthogonal  main-effect  plan  then  the  three-level-factor  identified  by 
column  and  the  three  two-level  factors  identified  by  columns  ^  ^ 

2  cannot  be  used. 


B.  Use  of  the  Catalogue 

In  this  section  wc  will  illustrate  the  use  of  the  catalogue  with  several 

examples. 

,10 

2_; 

The  index  of  orthogonal  main-effect  plans  given  in  secdon  D  of  this 
chapter  indicates  that  a  plan  can  be  obtained  for  th<*  2^^  experiment  in 
twelve  trials  from  basic  plan  4.  The  basic  plan  has  twelve  treatment 
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combinations  and  eleven  factors.  Choose  any  leu  columns  of  this  plan 
and  the  required  plan  is  obtained. 

(ii)  3x2^ : 

3 

The  plan  for  the  3x2  experiment  in  eight  trials  can  be  determined 
from  basic  plan  2.  The  footnote  to  this  plan  indicates  that  if  the  three- 
level  factor  is  cr.osen  then  the  two-level  factors  numbered  1,  2  and  3 
cannot  be  chosen.  Thus  the  plan  is  obtained  by  choosing  the  column 
representing  the  three-level  factor  and  any  three  of  the  four  columns 
4,  5,  6  and  7. 

3 

A  plan  for  the  3x2  experiment  in  nine  trials  is  given  by  basic 
plan  3.  The  plan  can  be  obtained  by  choosing  column  1  from  the  three- 
level  factors  and  columns  2,  3  and  4  from  the  two-level  factors.  It  is 
clear  that  a  plan  for  the  3x2"^  experiment  can  be  obtained  from  basic 
plan  3  by  choosing  any  one  of  the  four  columns  for  three-level  factors 
and  three  columr.3  from  the  two-level  factors,  no  column  of  the  two-level 
factor  having  the  same  column  number  as  the  column  number  of  the  chosen 
three-level  factv>r. 


(iii)  4^  X  3  X  2^ : 


The  index  indicates  that  an  orthogonal  main-effect  plan  for  the 
2  5 

4  x  3  X  2  experiment  can  be  constructed  in  sixteen  trials  from  basic 

plan  5.  The  plan  may  consist  of  the  columns  numbered  ^  and  ^  from 

3  111 

the  four-level  columns,  ^  from  the  three-level  columns  and  q, 

i  and  \  from  the  two-level  columns.  The  use  of  the  four-level 
3  4 

columns  1  and  f  eliminates  the  use  of  the  three-level  columns 
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i  2, 

identified  by  ^  and  ^  and  also  tha  two-level  columns  identified  by 

0  0  0  0  0,0,^^  -  *1.  .1.  Ill  ^  T  •  4. 

,  ^  and  y  .  The  use  oi  the  thrcie-level  column  ^  eliminates 

the  use  of  the  two-level  columns  2  2  and  2  - 

I  C  7 

(iv)  8^  X  4‘  X  : 

A  plan  for  the  8^x4^x2‘®  experiment  in  sixty-four  trials  may  be 
deduced  from  basic  plan  25.  If  the  three  eight-level  factors  chosen  are 
the  columns  identified  by  }.  ^  and  ^  then  the  four-level  factors 
identified  by  the  seven  columns  numbered 


0  0  0  0  0  0  0 

1  2  3  4  5  6  7 

$  !{c 

and  the  two-level  factors  identified,  by  the  columns  numbered  from  j  to 
2 

j  cannot  be  used.  We  then  can  choose  the  seven  four-level  factors  to  be 
the  columns  numbered 


0  1 

8  to  4  . 

*  * 

Thus  the  ten  two-level  factors  must  be  chosen  from  the  columns 
numbered  ^  to  ^  • 


C.  Tables  of  Orthogonal  Polynomials 

The  orthogonal  contrasts  which  define  effects  and  interactions  can  be 
readily  determined  from  a  table  of  orthogonal  polynomials.  The  ad¬ 
vantage  of  using  orthogonal  contrasts  to  define  effects  and  interactions 
arises  from  the  fact  that  orthogonal  polynomials  are  so  constructed  that 
any  term  of  the  polynomial  is  independent  of  any  other  term.  This 
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property  of  independence  permiis  one  to  compute  each  regression 
coefficient  independently  of  the  others  and  also  facilitates  testing  the 
significance  of  each  coefficient. 

Tables  of  Orthogonal  Polynomials  for  the  case  of  equally  spaced 
le'^els  are  readily  available,  e.g.  Fisher  and  Yates  (1938),  Anderson 
and  Houseman  (1942).  It  would  be  an  impossible  task  to  compute  the 
orthogonal  polynomials  for  unequally  spaced  levels.  However  a  oiinple 
procedure  for  computing  these  orthogonal  polynomials  is  available  and 
wiJl  be  presented  below.  If  equally  spaced  levels  do  not  each  occur  in  a 
plan  an  equal  number  of  times  the  published  tables  of  orthogonal 
polynomials  are  not  appropriate.  The  orthogonal  polynomials  for  equally 
spaced  levels  which  do  not  occur  in  a  plan  with  equal  frequency  must  be 
cornpxited  by  the  following  method  for  unequally  spaced  levels. 

For  any  set  of  orthogonal  polynomials  the  linear  contrast  is  of  the 

ft  -m  S(a  +  px)y^,  where  a  and  p  are  constants,  x  is  the  level  at 

which  the  factor  occurs,  is  the  response  to  the  treatment  combination 

with  the  factor  at  the  x  level  and  the  summation  is  over  every  value  of 

X  whicn  is  represented.  The  quadratic  and  cubic  contrasts  are  of  the 
2  2  3 

form  S(a  +  px  +  yx  )y^  and  S(a  +  px  +  yx  +  6x  )yj^>  respectively. 

The  extension  to  higher  order  contrasts  is  obvious.  Two  contrasts  are 
>  -thogonal  if  the  coefficients  of  each  contrast  sum  to  zero  and  the  sum  of 
products  of  the  corresponding  coefficients  of  the  two  contrasts  is  zero. 

We  will  illustrate  the  procedure  for  obtaining  orthogonal  polynomials 
for  unequally  spaced  levels  with  an  example. 

Consider  an  independent  variable  x  with  levels  0,  1,  2  and  4. 
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if  - 


I  ■, 


'i' 


The  coefficients  of  the  linear,  quadratic  a^d  cubic  ‘'ontTt  sts  for  this 
example  is  displayed  in  Table  16. 

TABLE  16 

COEFFICIENTS  OF  ORTHOJONAL  CONTRASTS 


Level  of 

X 

Linear 

Quadratic 

Cubic 

A 

V 

a 

a 

o. 

1 

G  T  p 

a  +  p  +  y 

r.  +  p  +  y  +  6 

1 

r  2 

a  +  2p 

a  +  2p  +  4y 

a  +  2p  +  4y  4-  86 

4 

a  -r  4p 

a  +  4p  +  I6y 

o  +  4p  +  I6y  +  646  j 

The  coefficients  of  the  linear  contrast  must  sum  to  zero.  Thus, 

4a  +  7p  =  0  . 

Setting  p  =  1  we  find  that  a  =  -  7/4.  In  order  that  the  coefficients  of 
the  orthogonal  contrasts  be  integers  reduced  to  lowest  terms  we  multiply 
these  coefficients  by  4  to  obtain  p  =  4  and  o  =  -  7 .  Substituting 
a  =  -  7  and  p  =  4  in  the  linear  contrast  given  in  Table  16,  gives  the 
linear  coefficients. 


Level  of 

X 


Coefficient  of 
lirear  contrast 
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The  coefficients  of  the  quadratic  contrast  must  sum  to  zero.  Hence, 


4a  +  7B  +  2iy  =  0  . 

The  sum  of  products  of  the  corresponding  coefficients  of  the  linear  and 
quadratic  contrasts  must  also  equal  zero.  Thus, 

35p  +  145y  =  0  . 


Solving  these  two  equations  to  obtain  integral  values  for  a,  P  and  y  we 
obtain  a  =  14,  p  =  -  29  a'd  y  =  7. 

Ji  we  substitute  these  values  in  the  quadratic  contrast  and  reduce  the 
resulting  coefficients  to  lowest  terms  the  coefficients  of  the  quad^•atic 
contrast  is  given  by 

Level  of  Coefficients  of 

X  Quadratic  contrast 


0 


7 


1 


-4 


2 


-8 


4 


5 


Sui.ilarly  the  sum  of  the  coefficients  of  the  cubic  contrast  and  the  sum  of 
products  of  the  corresponding  coefficients  of  the  linear  and  cubic  contrasts 
and  the  quadratic  and  cubic  contrasts  must  each  equal  zero.  Hence, 

4a  +  7p  +  21y  +  736  =  0 
35p  +  145y  +  5816  =  0 
44y  +  2526  =  0 

Solving  these  equations  to  obtain  integral  values  for  a,  p,  y  and  6  we 
obtain  a  =  -  36,  p  =  392,  y  =  -  315  and  6  =  55.  If  we  substitute  these 
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‘  •  > 


+  ▼  .  I/SI 


I-  kJ 


X  3‘ 


values  in  the  form  of  the  coefficients  of  the  cubic  contrast  given  in 
Table  16  and  reduce  tne  resulting  coefficients  to  lowest  terms,  the 
coefficients  of  the  cubic  contrast  are  given  by 


uevei  ci 


Coefficients  of 
Cubic  contrast 


The  orthogonal  polynomials  are  presented  in  the  following  table. 


TABLE  17 


ORTHOGONAL  POLYNOMIALS 


Level  of 


Linear 


Quadratic 


Cubic 


The  symbol  p  represents  one  unit  of  the  linear  effect  of  a  factor 
when  set  equal  to  unity.  In  order  to  obtain  integral  coefficients  p  was 
set  equal  to  4  and  hence  ^p  represents  one  unit  of  the  linear  effect. 
Consequently  the  linear  contrast  with  coefficients  given  in  Table  17 
represents  the  estimate  of  ^  the  linear  effect  of  the  factor.  It  is  easily 
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verified  that  the  coefficients  of  the  quadratic  contrast  are  given  by 


29 


X  + 


7 

2 


2 

X 


2 

where  x  =  0,  1,  2  and  4,  respectively.  Thus  the  symbol  jy 

represents  one  unit  of  the  quadratic  effect,  and  the  linear  contrast  with 

2 

coefficients  given  in  Table  17  represents  the  estimate  of  y  the  quadratic 
effect  of  the  factor.  Similarly  it  may  be  demonstrated  that  the  cubic 
contrast  wdth  coefficients  g’ ven  in  Table  17  represents  the  estimate  of 
12/55  the  cubic  effect  of  the  factor. 

This  constant  which  is  multiplying  each  effect  will  be  denoted  by 


and  in  the  tables  of  orthogonal  polynomials  the  value  of  A.  and  the  sum  of 
squares  of  the  coefficients  denoted  by  h,  will  both  be  given.  Thus  any 
contrast  defined  by  the  coefficients  given  in  the  tables  of  orthogonal 
polynomials  represents  ^  times  the  tippropriate  effect  of  the  factor.  It 
was  this  convention  by  which  the  coefficients  of  the  parameters  for  the 
example  in  section  H  of  Chapter  IV  were  calculated. 

In  the  tables  of  orthogonal  polynomials  the  coefficients  of  a  linear 
contrast  will  be  denoted  by  Oj,  the  coefficients  of  a  quadratic  contrast 
by  0^  and  so  on.  The  levels  of  a  factor  will  be  denoted  by  x. 


112 


170  1,326  78  {  476  1,972  58  828  4^692 


116 


11,424  12  14,304 

180  1,120  125,970  72  684  288  527,460 


1 
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168  168  264  (  4,284  58,548 


2,772  990 


D.  Index  of  Orthogonal  Main-Effect  Plans 


The  index  presented  in  tliis  section  indicates  the  basic  plan  from 

which  any  orthogonal  main- effect  plan  can  be  deduced  with  a  minimum 

number  of  trials.  Unless  the  experimenter  must  use  a  plan  with  a 

minimum  number  of  trials  there  is  usually  a  choice  of  basic  plans  from 

which  an  orthogonal  main-effect  plan  can  be  constructed.  For  example, 

4  3 

the  basic  plan  from  which  the  main-effect  plan  for  the  3  .  2  experiment 

can  be  constructed  in  16  trials  is  basic  plan  5.  However  a  plan  for  this 

experiment  can  be  constructed  from  basic  plan  7  in  18  trials.  The  use 

of  one  basic  plan  over  another  depends  on  which  contrasts  are  deemed  to 

4  3 

be  most  important.  The  orthogonal  main- effect  plan  for  the  3  .  2 
experiment  in  l6  trials  estimates  the  two-level  factor  with  an  efficiency 
of  unity  and  estimates  the  linear  effect  of  the  three-level  factors  with 
efficiency  3/4  and  their  quadratic  effects  with  efficiency  9/8.  The  plan 
with  18  trials  estimates  the  effects  of  the  three-level  factors  with  an 
efficiency  of  unity  and  the  effect  of  each  two-level  factor  with  an  efficiency 
of  8/9.  If  the  effects  of  the  three -level  factors  are  the  more  important 
then  the  plan  with  18  trials  should  be  chosen  and  if  the  effects  of  the  two- 
level  factors  are  the  more  important  then  the  plan  with  16  trials  should 
be  chosen. 

The  notation  used  in  the  index  requires  some  explanation.  The  plan 
20+n, 

3  *.2  ,  En.  =  5  in  54  trials  indicates  that  orthogonal  main-effect 

plans  can  be  constructed  in  54  trials,  from  basic  plan  22,  for  the 
following  experiments: 
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3^^2^,  3^^2^  3^12  and  3^^. 


The  notation  is  used  to  reduce  the  number  of  entries  necessary  to  list  all 

the  possible  plans.  The  plan  4^.  3*^.  n  =  0, 1,  2  in  32  trials 

indicates  that  the  plans  for  the  4^.  2^^  ,  4*^.  3.  2}  and  4"^.  3^.  2^ 

experiments  can  be  constructed  in  32  trials.  The  plan 
t,  t^  n,  n, 

6'^5  .4  .3“. 2  ,  St.  =  8,  Sn.  =  8-4nj  represents  the  following 


experiments: 
h  ‘2  ,^^2^3 


tj^  t^  n^  n.. 


6  S  ^.3  ^2  St.  =  8,  Sn.  =  8;  6  ^  5  4.  3  ^2  St.  =  8,  Sn.  =  4; 


^1  ^2  2 
6^5.4'^,  St.  =  8. 

i 
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Plan  Number  of  Basic  Page 

trials  plan 


Z^ 

4 

1 

139 

2'^ 

8 

2 

139 

,11 

12 

4 

140 

^15 

16 

5 

141 

IQ 

2"  ' 

20 

8 

142 

24 

9 

143 

2^^ 

28 

12 

146 

2^' 

32 

13 

147 

2^^ 

36 

15 

149 

2^9 

40 

16 

150 

2"^ 

44 

17 

151 

2"^ 

48 

18 

152 

2^' 

52 

21 

155 

2^5 

56 

23 

157 

2^9 

60 

24 

158 

2^3 

64 

25 

159 

3.2^ 

8 

2 

139 

3.2^“ 

16 

5 

141 

3.22» 

32 

13 

147 

3.  ,60 

64 

25 

159 

3^.2" 

9 

3 

140 

3^.29 

16 

5 

141 

IZZ 


I 


Plan 

Nvimber  of 
trial 9 

Basic 

plan 

Page 

32.  2“ 

27 

11 

145 

32  2-5 

32 

13 

147 

32,257 

64 

25 

159 

2 

9 

3 

140 

35.  2^ 

16 

5 

141 

35.  2^° 

27 

1 1 

145 

35. 2^2 

32 

13 

147 

,3  ,54 
j  .  z 

64 

>5 

159 

3^ 

9 

3 

140 

3^  2/ 

16 

5 

141 

3^.  2^ 

27 

11 

145 

3^  2^9 

32 

13 

147 

3^  2^^ 

54 

22 

156 

3^2^^ 

64 

25 

159 

35 

16 

5 

141 

35.22 

18 

7 

142 

35.28 

27 

11 

14S 

35.2^6 

32 

13 

147 

35. 220 

54 

22 

156 

35.  2^8 

64 

25 

159 
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Plan 


Number  of 
trials 


Basic 

plan 


Page 


3^.  2 

18 

7 

142 

3^2'^ 

27 

11 

145 

3^2^^ 

32 

13 

147 

3^  2’-^ 

54 

22 

156 

3^  2^5 

64 

25 

159 

3^ 

18 

7 

142 

3^.  2^ 

27 

11 

145 

3^2^° 

32 

13 

147 

3'.2^° 

54 

22 

156 

3^. 2^2 

64 

25 

159 

3®.  2® 

27 

11 

145 

3®.  2“^ 

32 

13 

147 

3®.  2^^ 

54 

22 

156 

3».2” 

64 

25 

159 

3^.  2^ 

27 

11 

145 

3^.2^® 

54 

22 

156 

3’.  3’‘' 

64 

2.5 

159 

3*".  2^ 

27 

11 

145 

3'®.2'5 

54 

22 

156 

s'". 2^^ 

64 

25 

159 
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Plan 


Number  of 
trials 


Basic 

plan 


Page 


3^^.  2^ 

27 

11 

145 

2^^ 

54 

22 

156 

311^30 

64 

25 

159 

3^^  2 

27 

11 

145 

54 

22 

156 

3^1  22-^ 

64 

25 

159 

3I2  ^28 

81 

26 

162 

3*3 

27 

11 

145 

3^3.212 

54 

22 

156 

3I3  224 

64 

25 

159 

3^3.227 

81 

26 

162 

3^1  2^^ 

54 

22 

156 

3I4  ^21 

64 

25  . 

159 

3I4  226 

81 

26 

162 

3^3.2'° 

54 

22 

156 

3I5  218 

64 

25 

159 

3^5.  2-^3 

81 

26 

162 

336.  2^ 

C  A 
»/”X 

22 

156 

3I6  2I5 

64 

25 
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BASIC  PLAN  14:  8  ;  7  ;  6  ;  5  ;  4^  ;  3^  ;  2^'^  ;  32  trials 
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